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1.0  SUMMARY 


The  overall  preliminary  theoretical  flutter  analysis  of  the  XV-5A 
Lift  Fan  Aircraft  is  presented  in  three  volumes.  This  is  Volume  I 
which  reports  the  investigation  of  the  flutter  characteristics  of 
the  wing.  Volumes  II  and  III  detail  the  analysis  of  the  empennage 
and  control  surfaces  respectively. 


The  wing  theoretical  flutter  analysis  was  a  joint  effort  by  Ryan 
Aerospace  and  Computer  Engineering  Associated  (CEA),  an  affiliate 
of  Susquehanna  Sciences,  Inc. ,  Pasadena,  California.  The  program 
occupied  the  CEA  direct  analog  computer  during  the  period  of  3  May 
to  25  May,  1962.  The  primary  concern  of  the  program  was  that  of 
flutter  evaluation  of  the  XV- 5A  wing  in  the  conventional  flight 
mode.  In  addition,  vibration  mode  shapes  and  frequencies  were 
obtained  for  several  wing  configurations  (stiffness  variations) 
both  with  and  without  gyroscopic  forces  due  to  the  rotating  fans 
in  the  wings. 


Results  of  the  overall  CEA  analysis  have  indicated  that  the  XV-5A 
aircraft  (wing),  on  the  basis  of  the  preliminary  data  utilized  in 
the  analysis,  is  free  of  flutter  within  the  proposed  flight 
envelope.  Studies  have  shown  that  mass-balancing  of  the  aileron  is 
a  critical  factor,  and  must  be  carefully  evaluated  for  future  design 
efforts  on  the  XV-5A  aircraft.  Similarly,  the  overall  spring 
restraint  (Kf}}  does  affect  the  flutter  characteristics,  but  the 
variation,  at  most,  will  still  lie  above  15%  on  V^  or  M^.  The 
flutter  mode,  that  of  antisymmetric  wing  bending  coupled  with 
aileron  rotation,  indicates  that  proper  aircraft  simulation  must 
be  used  in  future  flutter  evaluation  of  the  XV-5A  aircraft  with 
due  attention  given  to  fuselage  and/or  aircraft  degrees  of  freedom. 
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BLANK  PAGE 


2.0  INTRODUCTION 


This  report  is  the  final  presentation  of  the  preliminary  wing  flutter 
and  vibration  (modal)  analysis  of  the  U.S.  Army  XV-5A  Lift  Fan 
Research  Aircraft,  illustrated  in  Figure  1.0  performed  on  the  CEA 
direct  analog  computer.  The  XV- 5A  is  a  V/STOL  aircraft  designed  for 
research  flight  testing  of  the  General  Electric  X353-5  Lift  Fan 
Propulsion  System. 

The  use  of  a  direct  analog  computer  such  as  the  CEA  computer  was  dic¬ 
tated  by  close  scheduling  called  for  on  the  XV- 5A  Program  and,  conse¬ 
quently,  required  quicker  answers  relating  to  flutter  evaluation  than 
would  have  been  possible  via  conventional  means,  (e.g.  solution  of  the 
appropriate  equations  of  motion  via  digital  computer  techniques).  As 
stated  in  Reference  2,  a  direct  analog  computer  offers  an  advantage 
over  other  computational  procedures,  since  it  permits  the  analyst  to 
observe  immediately  the  effect  of  design  changes  on  the  solution. 

The  formulation  of  an  aeroelastic  problem  for  solution  by  the  direct 
analog  computer  is  accomplished  by  first  representing  electrically 
the  idealized  mechanical  model  of  the  structure  under  investigation, 
i.e.  the  circuit  is  a  direct  analogy  of  the  structure.  Next,  the 
external  forces  acting  on  the  structure,  (in  this  case  aerodynamic 
forces),  are  simulated  by  an  electronic  circuit  which  generates  an 
electrical  analogy  for  the  aerodynamic  forces  as  a  function  of 
structural  deflections  and  solution  parameters,  such  as  aircraft 
velocity.  Combining  both  the  circuit  which  is  analogous  to  the 
structure  to  that  which  is  analogous  to  the  aerodynamic  forces, 
results  in  an  analog  simulation  of  an  airborne  elastic  vehicle. 
References  3  through  5  provide  additional  references  related  to  the 
use  of  direct  analog  computers. 

Equations,  following  the  technique  of  Reference  6,  are  developed  in 
Section  3.0,  Method  of  Approach.  These  equations  are  presented  as 
an  analogy  to  the  actual  flutter  and  vibration  analysis  performed 
on  the  CEA  direct  analog  computer. 

It  is  to  be  emphasized  that  these  equations  are  not  solved  by  the 
direct  analog  computer,  and  are  offered  only  as  a  mathematical  guide 
for  an  understanding  of  the  analysis. 
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3.0  METHOD  OF  APPROACH 


Of  the  several  methods  of  determining  the  flutter  characteristics  of  lift¬ 
ing  surfaces,  the  influence-coefficient  method  of  analysis  for  calculating 
the  response  of  a  flexible  surface  to  an  oscillating  disturbing  force  is 
similar  to  that  performed  by  direct  electronic  analog  treatment.  The 
commonly  used  Rayleigh-Ritz  technique,  (formulating  the  equations  di¬ 
rectly  in  terms  of  chosen  modal  functions),  is  essentially  redundant  for 
the  method  of  approach  chosen  for  flutter  analysis  of  the  wing. 

The  deflection  q  of  a  point  x,y  on  a  structure  subjected  to  r.  set  of  dis¬ 
crete  forces  or  moments  acting  normal  to  the  surface  (assumed  lying  in 
the  x-y  plane  and  acting  in  the  direction  of  z  may  be  expressed  as 

{q}  =  [ C  J  {Q}  (1) 

Equation  1  is  the  matrix  form  of  a  set  of  simultaneous  equations  where  an 
element  Cjj  relates  the  generalized  displacement  q.  to  the  generalized 
force  Qj.  In  structural  nomenclature,  the  element  c^j  is  a  flexibility 
influence  coefficient  which  expresses  the  linear  or  angular  deflection  at 
point  i  due  to  the  unit  force  or  moment  acting  at  point  j.  The  use  of 
generalized  displacements  and  forces  in  Equation  1  implies  that  these  may 
be  either  linear  and/or  angular  deflections  and  corresponding  forces 
and/or  moments  with  appropriate  flexibility  influence  coefficients. 

In  general,  the  forces  to  be  considered  in  Equation  1  are  an  impressed 
sinusoidal  force,  inertia  forces,  structural  damping  forces,  aerodynamic 
forces  and  gyroscopic  forces. 

Consider  now  each  of  the  above  forces: 


Impressed  Sinusoidal  Force 


The  impressed  sinusoidal  force  considered  herein  is  simply  a  discrete 
force  or  moment  applied  at  a  node  (point  on  the  surface  in  the  x-y  plane) 
used  for  purposes  of  determining  the  flutter  and  modal  response,  as  will 
be  discussed  later.  In  general,  the  force  may  be  expressed  as 


Q,  =  Q  .  e 
j  oj 


icot 


(2) 


where  the  subi  cript  j  denotes  the  node  at  which  the  force  or  moment  is 
applied. 
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Inertia  Forces 


The  inertia  forces  considered  herein  are  the  forces  set  up  by  consequence 
of  the  accelerations  developed  in  the  structure.  Since  for  flutter  and 
modal  analyses,  we  are  concerned  with  harmonic  motion,  the  inertia 
forces  may  be  expressed  as 

{Q}  Jj{q}  =U2e‘"‘  rjJ{qo}  (3) 

where  o>  is  the  angular  frequency  and  the  matrix  [  J  ]  is  .1  diagonal 
matrix  of  masses  and  mass  inertias  lumped  at  the  nodes. 

Structural  Damping  Forces 

The  structural  damping  forces  acting  in  the  structure  are  forces  in  phase 
with  the  velocity  and  proportional  to  the  restoring  force.  Therefore,  we 
may  write  for  the  damping  forces, 

{Q>  =  -g/w  [K]{q}  =  -igela*  [Kj{qo}  (4) 

where  g  is  the  coefficient  of  structural  damping  and  the  matrix  I  K  J  is 
the  matrix  of  stiffness  influence  coefficients.  The  stiffness  matrix  shown 
in  Equation  4  is  the  inverse  of  the  flexibility  matrix  shown  in  Equation  1, 
i.e.,  [  K  J  =  l  C]'1. 

Aerodynamic  Forces 

The  aerodynamic  forces  to  be  considered  herein  may  be  represented  by 
an  expression  also  called  an  influence  coefficient,  although  it  is  of  a 
different  nature  than  that  denoted  by  Equation  1.  Therefore, 

{Q  }  =  [  A  J  {q }  =  eiWt  (A]  {qQ}  (5) 

where  the  matrix  [  A  ]  contains  aerodynamic  influence  coefficients  whose 
element  A^  represents  the  aerodynamic  force  and  moment  induced  at 
point  i  by  a  unit  deflection  at  point  j.  The  above  matrix  [A  J  is  a  complex 
matrix  having  real  and  imaginary  parts. 


Gyroscopic  Foices 


The  gyroscopic  forces  (moments)  developed  on  the  wing  structure  are  due 
to  the  wing  fan  and  may  be  expressed  as 


[ a  j{qk } 


icot 

iu>e 


[G] 


(6) 


where  the  subscript  k  denotes  the  node  at  which  the  gyroscopic  moments 
act  and  the  matrix  [  G  J  is  simply  a  matrix  expressing  the  angular 
momentum  of  the  rotating  mass. 


Substituting  Equations  2  through  6  into  Equation  1  results  in  the  following 
matrix  equation  for  forced  vibration  in  an  airstream: 


{q  J  =  |C] 


Q 


2 


+  W  [C]h  JvJiq^}  -  ,K{qo>  +  I  C  ]  I  A  J  {  qQ} 


+  ,W[C](G] 


(7) 


Since  Q  =  qQ.,  where  J.  is  either  the  mass  or  mass-inertia  lumped 

°j  J  j  J 

at  point  j,  Equation  7  may  be  written  as 


2  2 
{q  }  ="  l  C  )t~  JJiq}  -  iB{q  }  +"  ICir  J  J 


o 


o 


+  [c  ]  [  A  ]{qQ} 


+  iwIcHg] 


(8) 


where  the  column  matrix  Jq0^ 


is  composed  of  zeros  plus  one  finite 


element  which  is  q^,  the  displacement  at  point  j  due  to  the  impressed 

sinusoidal  force  Q0^  and  the  column  matrix  |  qG^j  is  composed  of  zeros, 

plus  two  finite  elements  which  are  the  angular  displacements  about  the 
x  and  y  axes  at  point  k. 
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Rearranging  Equation  8  results  in  an  expression  for  the  column  of 
generalized  displacements: 


W 


2  2 
"  ICJCVU  -  l*UI  +  ICHA)  {qj  +"  (CjhJj 


+  ‘“iCl  IG) 


2  2 
II)  -  "  [C]  fJj  -  (Cl  (A]  +  **11]  {qj  -  "  [Cl  M 


V 


+  l“lc]  IG) 


V 


Dividing  through  by  of 


^(Il-ICl  W-^(C1  (A)  +  III 

or  a)  a) 


+  [Cl  (Gj 

Cl) 


=  1C]  MJ 


% 


or,  letting  ft  =  \  (1  +  ig)  and  [Aj  =  o  [A] 


CO 


“(I)  -  (Cl  Ml  -  "(C11A) 


{qj  -  (ClfJj 


V 


+  7:10  (Gl 
CO 


V 


and 


V  ■ 


o 


[IJ-  (Cl  fJj  -  a[C]  [A] 


-1 


icj  rjj 


+  ~IC]  [GJ 

CO 


Equation  9  is  written  in  a  general  sense,  and  is  applicable  to  either 
modal  analysis  or  flutter  analysis,  as  will  be  shown  in  subsequent 
sections. 
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3.1  IDEALIZATION  OF  AIRCRAFT 

The  technique  of  performing  both  a  flutter  and  vibration  analysis  of  an 
aircraft  usually  requires  an  idealization  of  the  basic  structure  and  its 
external  loading  which  is  amenable  to  analytical  treatment.  The  ideali¬ 
zation  must  be  of  sufficient  rigor  to  meet  the  demands  of  engineering 
accuracy,  yet  be  consistent  with  the  capability  of  available  digital  or 
analog  computers.  Reference  1  presents  the  idealization  as  used  in  the 
Computer  Engineering  Associates'  flutter  and  vibration  simulation. 

However,  for  purposes  of  completeness,  the  discussion  will  be  repeated 
here  with  added  commentaries  on  the  applicability  of  passive  analog 
computers  to  the  present  problem. 

3.1.1  Structural  Representation 

Since  the  wing  of  the  XV-5A  aircraft  is  of  low  aspect-ratio  (AR=3.4)  and 
of  unusual  structural  arrangement,  (i.e. ,  with  a  large  cutout  for  wing 
fan),  it  was  decided  to  consider  the  wing  more  as  a  deforcing  plate  with 
a  possible  elastic  curvature  of  appreciable  magnitude  in  two  directions, 
(spanwise  and  chord-wise).  Secondly,  the  root  restraint,  that  of  two 
spars  with  a  forward  nose  box,  presented  peculiarities  which  play  an 
important  part  in  the  deformations  of  the  surface.  For  these  reasons,  it 
was  also  decided  to  consider  the  wing  as  an  assemblage  of  a  finite  num¬ 
ber  of  elastic  components,  e.g.,  flexural  members,  torque  tubes  and 

torque  boxes.  The  analysis  was  primarily  concerned  with  the  wing  dynamic 
characteristics,  with  allowance  for  fuselage  and  empennage  effects.  How¬ 
ever,  these  latter  structures  were  treated  as  simple  beams,  subject  to 
engineering  beam  theory,  in  order  to  keep  within  the  capabilities  of  the 
analog. 

The  following  discussion  interprets  each  component.  These  were  then 
set  up  on  the  analog  computer  by  passive  elements. 

3. 1.1.1  Beam  Elements  (Flexure  Members) 

The  beam  elements  of  the  idealized  structure  can  be  divided  into  four 
categories.  These  are: 

1.  Beams  flexible  in  bending  and  shear 

2.  Beams  flexible  in  bending  but  rigid  in  shear 

3.  Beams  flexible  in  bending  but  with  no  shear  carrying 

capability  (no  shear  web) 

4.  Rigid  beams 
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Beams  of  Types  2,  3  and  4  are  special  cases  of  flexibility  of  the  general 
beam  element  of  Type  1 . 

Consider  now  the  beam  shown  in  Figure  2  which  represents  the  beam 
element  of  Type  1  above.  The  beam  is  deformed  by  forces  and  moments 
applied  at  panel  points  1  and  2.  Panel  points  1  and  2  represent  the  nodes 
of  a  lumped  parameter  system  resulting  from  the  idealization  of  the 
actual  continuous  structure. 

The  load-deflection  relationships  of  the  beam  may  be  expressed  by  the 
following  matrix  equation, 


e„ 


M. 


2f  1  (l-*//)2  r 

1  I 


o 


dx 

Gth 


1 

l 

f  dx 

1 

1  0  1 
_l_  _ 

J  El  dX 

0 

1 

- 1- 

J  El 

0 

-1 

_l_ 

0 

|  0  0 
““  1  ■■  ■ 1  “ 

/  d-x//) 
tj  — rr  dx 

o 


El 


-1 


1  i 


0  0 


M 


2 


0 


(10) 


The  above  equation  was  evaluated  with  the  assumption  that  the  flexibility 
(l/EI)  varies  linearly  over  each  beam  segment  and  that  the  beam  height 
(h)  was  constant  over  the  segment. 


3 . 1 . 1 . 2  Torque  Tube  Elements 

The  torque  tube  elements  as  shown  in  Figure  3  of  the  idealized  structure 
may  be  expressed  by  the  following  matrix  equation, 


dx 

GJ 


1 


I 

I 

T 

I 

I 


1 

f  1 

Mn 

2 

0 

0, 

1 

k  4 

(11) 
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The  equation  was  evaluated  with  the  assumpti  a  that  the  flexibility 
(l/GJ)  varies  linearly  over  each  torque  tube  element. 

3. 1.1.3  Torque  Box  Elements 

The  torque  box  is  best  described  as  an  element  composed  of  two  shear 
panels,  symmetrically  located  above  and  below  the  neutral  plane  of  the 
structure.  Assuming  that  loads  and  displacements  of  the  shear  panels 
are  antisymmetric  about  the  neutral  plane,  then  shear  loads  on  the 
boundaries  of  the  panels  are  moments  about  the  neutral  plane  and  dis¬ 
placements  of  the  panels  are  rotations  about  the  neutral  plane. 

Considering  a  typical  shear  panel  as  shown  in  Figure  4,  the  basic  equa¬ 
tions  are: 


For  equilibrium, 


i  cos  6  1  cos  6  ’  d  d„ 

1  2  4  2 


(12) 


The  strain  energy  stored  in  the  panel  is, 


12  l 

U  =  -  F  - 

2  4  Gtd, 


+  tan  <52) 


2d. 


l+  aaf^(tan"6i+  tan*i  tan62 


(13) 


3.1.2  Aerodynamic  Representation 


For  analog  computation,  simulation  of  the  unsteady  aerodynamics  for  a 
flutter  analysis  of  a  lifting  surface  requires  the  use  of  assumptions  or 
empirical  approximations.  As  stated  in  Reference  1,  reasonably  relia¬ 
ble  methods  are  available  for  the  calculation  of  steady  state  aerodynamic 
loads.  However,  due  to  the  greater  complexity  of  the  nonstationary 
problem,  available  methods  for  computing  oscillating  aerodynamic  forces 
are  less  exact.  Therefore,  the  aerodynamic  representation  used  in  the 
flutter  analysis  of  the  XV-5A  aircraft  was  accomplished  by  a  rational 
combination  of  the  following: 


A.  Mathematical  expressions  derived  for  incompressible  flow  over 
an  infinite  aspect  ratio  surface. 


n 


B.  Steady  state  aerodynamic  influence  coefficients  relating  force 


and/or  moments  to  linear  and/or  angular  deflections,  calculated 
on  a  digital  computer. 

C.  Various  approximations,  proved  reasonable  in  the  general  flutter 

analysis  experience. 

3. 1.2.1  Development  -  Wing 

With  reference  to  statement  "A"  of  Paragraph  3. 1.  2 
motion  of  a  rigid  chord  is  defined  (see  Figure  5): 

the  following  wing 

• 

a  =  a  +  h/V 

(14) 

q  =  c/V  at 

(15) 

qp  =  c/V  P 

(16) 

where: 

V  =  Airstream  velocity 

c  =  Aerodynamic  chord 

Of  =  Angle  of  rotation  (+  L.  E.  up) 

h  =  Vertical  airfoil  velocity  at  a  reference  axis  (+  down) 
a  =  Angle  of  attack  (at  the  specific  reference  axis) 
q  =  Dimensionless  pitch  velocity 

/3  =  Angle  of  rotation  of  control  surface  (+  T.E.  down) 

x  =  Distance  from  leading  edge  to  moment  axis  as  a  fraction  of  the 

total  chord  (c) 

x^  =  Distance  from  leading  edge  to  rotation  center  as  a  fraction  of  the 
total  chord  (c) 

qp  =  Dimensionless  control  surface  pitch  velocity 

x  =  Control  surface  hinge  location  with  respect  to  the  leading  edge 
3- 


*m  i  * 


From  Equations  72,  73  and  74  of  Reference  7,  the  lift,  moment  and  con¬ 
trol  surface  hinge  moment  per  unit  span  may  be  written  as 


(t'jv2)<c> 

C  *  a 
L 

t  +  , 

CL  *“ 

+  , 

C  *  13 

\  / 

a 

q 

P 

M 


P 


-1 

C  *  a 

Ma 

>  +  < 

I 

C  *  q 

M  H 

q  J 

’  +  * 

p 

M 

P  I 

(17) 


where  the  asterisk  between  the  coefficient (C,  ,  etc.)  and  the  particular 

'  i O' 

coordinate  (a,  etc.)  represents  the  indicial  response,  i.e.,  the  coeffi¬ 
cients  (Ct  ,  etc.)  represent  the  indicial  response  to  a  unit  step  change 

in  the  particular  coordinate  involved.  The  indicial  response  of  a  physi¬ 
cal  system  is  the  time  variation  of  response  to  a  step  input  if  the  system 
is  initially  at  rest. 


Following  the  procedure  developed  in  Reference  7,  Equations  17  through 
19  may  be  written  as 


L 

=  L  +  L  +  L  +  L 

12  3  4 

(20) 

M 

=  M  +  M_  +  M„  +  M  4 

12  3  4 

(21) 

R 

■  "l  ‘  "2  +  »3  +  »4 

(22) 

where  the  subscript  1  stands  for  the  circulation  components;  2  for  the 
camber  components;  3  for  the  first  derivative  impulsive  terms  and  4  for 
the  equivalent  inertia  terms. 
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Therefore,  as  in  Reference  7,  the  components  are: 
Circulation  Components  (Xq  -  1/4;  =  3/4) 

Ll  =  ~(tfiv2)  (C)  <Zn)  |C(k)  ("  +  V  +  V+  W2  I  4| 
Mx  =  0 

H1  =  W3  <C>  L1 


Camber  Components 


l2  =  0 


m2  = 


»2  ■ 


1  „2'lc 

2Py 


y(i) 


1  v2'lc 
2ov 


a  +  W  -/3+W  8 
4  c  5 


we"  +W7  o*+V 


First  Derivative  Impulsive  Terms 


m3  = 


H3  = 


)-Xo)“+{WU  +  (Xa-V  V* 


Equivalent  Inertia  Terms 

2 


L4  = 


7T  P  C 


h  +(  g  -X1)c"  +  w10c  0 


(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 
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(33) 


where  the  W  functions  are  all  functions  of  X  and  are  similar  to  Theo- 

ji 

dorsen's  T  functions,  but  expressed  as  explicit  functions  of  the  parameter 
<p  =  cos-1  (2Xa-l). 


With  reference  to  Statement  "C"  of  Paragraph  3. 1.2,  the  following 
approximations  will  be  made  to  the  control  surface  "W"  coefficients: 


W2/  =  «W./  =  4  Wo/ 

2/w  *  5/w.  *  8'w_ 

14  7 


W  6  =  ' 2  W5 

Numerical  values  of  the  above  constants  are  presented  in  Reference  7, 
and  the  above  approximations  are  reasonable  for  the  range  of  values  of 
Xa  to  be  used  in  the  analysis.  Table  1  lists  the  possible  range  of  Xa  with 
the  appropriate  W  relationships. 

The  following  effective  control  surface  rotation  is  now  defined: 

W 

0  =  0  +  %  2/Wx  0  (35) 


Substituting  the  above  approximations  and  Equations  35  and  14  into  the 
circulation  and  camber  components,  (i.e.,  Equations  23  through  28), 


C(k)  (a  +  Wx  0 ) 


0 


(36) 

(37) 


H1  =  W3  <C)  L1 


(38) 


(39) 
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(40) 


M2  -  -(lpv2)<c2>  (I)1V  “+ W4  °' 

“2  -  -(4pv2)(c2>  (l)t^  ° 


2  W  —  4 
5  V 


.2  Wg  «  +  W?  /?  +  0.6  Wg-/3 


The  total  circulatory  and  camber  terms  then,  are: 


L  =  L  +  L 
12 

=  -(±pv2\  (C)  (2ir){c(k)  [5  +  Wl  ?]} 


M  =  M  +  M 


=-(ipv2)<c2)(f)l7° 


+  W  B  +0.2  W  -0 
4  5  V 


H  = 


VH2 


-(i '**)'■■>  (i)I 


(41) 


(42) 


(43) 


16W  C(k)  a  +  (16W  C(k)  W,  +  W_)  £ 
O  O  17 


+  0.2  W  7,  <*  +  0.6WQ  £/? 

D  V  o  V 


(44) 


The  first  derivative  impulsive  terms  (Equations  29  through  31)  are  known 
to  "smear"  out  due  to  compressibility  effects.  At  the  speeds  to  be  con¬ 
sidered  in  the  analysis,  these  terms  were  neglected. 

The  equivalent  inertia  terms  (Equations  2  through  34)  define  "aero¬ 
dynamic  mass"  which  is  added  to  the  structural  mass  of  the  surface. 
Since  it  is  small  relative  to  the  mass  of  the  wing,  these  terms  are 
neglected. 

Now  with  reference  to  Statement  "B"  of  Paragraph  3.1.2,  aerodynamic 
influence  coefficient  matrices  (QJ,  will  be  defined  in  steady-state  level 
flight  as  follows: 


For  the  lift  at  the  i^  section  (strip)  in  pounds, 
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For  the  moment  at  the  i  section  about  the  aerodynamic  center  in  inch- 
pounds, 


{M(}  =-(ipV2)cM^  IQ'fiOij}  (46, 

For  the  control  surface  hinge  moment  in  inch-pounds, 


In  the  above  expressions,  the  subscripts  i  and  j  refer  to  the  strips  into 
which  the  wing  has  been  divided  for  purposes  of  aerodynamic  simulation. 

Equations  45  through  47  define  the  steady-state  aerodynamic  forces  on 
the  surface  and  are  presumed  to  include  the  effects  of  spanwise  coupling 
between  strips  (i.  e.,  influence  of  strip  j  on  strip  i),  compressibility, 
aerodynamic  balance  on  the  aileron,  etc. 

Looking  now  at  the  development  of  Equations  45  through  47,  the  aero¬ 
dynamic  matrices  etc.  may  be  defined  as  weighted  numbers 

which  when  multiplied  by  the  value  of  the  effective  angle  of  attack,  etc . 
at  the  j**1  station,  will  give  the  contribution  of  the  j**1  station  to  the  i^ 
station  angle  of  attack,  etc.  Integration  of  the  running  lift,  etc.  is 
implied  in  the  above  influence  coefficient  matrices.  However,  for  pur¬ 
poses  of  the  following  discussion,  an  aerodynamic  influence  coefficient 
matrix  will  be  defined  as  a  matrix  which  relates  local  lift  coefficient  and 
local  angle  of  attack.  These  coefficients  may  be  thought  of  as  operators 
which  have  the  property  of  being  able  to  operate  on  an  arbitrary  angle  of 
attack  distribution  to  produce  the  resultant  span  loading. 

The  steady-state  loads  for  the  analysis  were  derived  on  the  basis  of 
strip  theory  and  also  by  the  lifting-line  concept  of  treating  aerodynamic 
loading  (i.e.,  the  distributed  loading  is  replaced  by  a  spanwise  lift,  and 
moment  distribution  assumed  to  act  along  the  line  of  chordwise  centers 
of  pressure). 

The  lifting-line  aerodynamic  influence  coefficient  matrix  may  be  defined 
as  follows: 

=  IQ)  {a,  fi}  (48) 
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where  is  the  local  sectional  lift  coefficient  and  c  is  the  local  chord. 
The  term  cc 4  /c  represents  the  nondimensional  loading  resulting  from 
an  arbitrary  angular  deflection  ot  or  p  distribution. 

Consider  now  the  following  loading  due  to  an  angle  of  attack  distribution 
and  due  to  a  trailing  edge  flap  rotation  distribution: 


The  loading  due  to  an  angle  of  attack  is: 


= 

y 

+ 

Yf 

Wa 

*ot 

(49) 


where  each  individual  loading  may  be  expressed  in  terms  of  influence 
coefficients  as, 


w 


a 


“  iV1 


(yJ  =  ‘V  {a} 


(50) 


(51) 


The  loading  due  to  a  flap  rotation  is: 

<y  ■  {v+{v  <52> 

and  expressing  each  component  in  terms  of  influence  coefficients  as, 

{W=,V{'J}  <53) 

{  y  -  iQgj)  {*}  (54) 

The  total  loading  on  the  wing  is  therefore: 

{y  }  =  { v  }  +  {y  } 

w  wot  wj3 

=  |vw+|vw  (55) 
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The  running  lift  in  pounds  per  unit  span  may  be  written  as: 


{fw}  =  -cq{Vw}  (56) 

The  total  lift  force  acting  at  each  of  the  specified  panel  points  is: 

{hi  =  r  R  }  (57) 

w  w 

where  f~  H  -J  is  simply  an  integrating  matrix, 
w 


Substituting  Equations  56  and  55  into  Equation  57  results  in  a  final 
expression  for  the  lift, 

{L}  =  -cqr-Kv_]  lQw„Haj-cq  rHwJ[Qw/iJ{/)}  (58) 

The  (Q  ]  and  [Q  n]  matrices  were  computed  by  the  method  of 
„  WfV  ,  wp 
References  8  and  9. 


The  loading  on  the  flap  due  to  an  angle  of  attack  and  to  a  flap  deflection 
was  computed  on  a  two-dimensional  (strip  theory)  basis.  Accordingly, 
we  may  write  for  the  resultant  hinge  moment  in  Ibs-ft.  per  unit  span, 


h 


a/3 


0  + 


ac. 

_ h 

3  a 


a 


(59) 


For  all  stations, 

{hi  =  -q  retire,  .]{<*}  -q  r-cf  -J('C  (60) 

a  hq  a  h0 

17ie  total  moment  acting  at  each  of  the  specified  panel  points  is: 


The  loading  on  the  wing  due  to  a  flap  deflection  also  gives  rise  to  a 
moment  about  the  panel  aerodynamic  center  (a.c.).  Accordingly,  by 
analogy  with  Equation  48,  we  may  write  the  loading  (moment)  produced 
by  a  flap  rotation  as, 


2  „ 
c  C 


{y}  = 


■  ‘V  <*> 


Therefore,  the  running  moment  about  the  a.c.  is: 
{m}  =  -  q  c  { y}  =  -  qc  [ Q^]  {#} 


And  the  total  moment  acting  at  each  panel  point  is: 

{m}  =  -qc  rRwJ  IQW/J]  {/»}  (64) 

Computation  of  the  matrix  product  T  R^--]  lQw^l  was  done  by  the  method  of 
Reference  9,  and  therefore  was  based  upon  the  lifting-line  concept. 

Comparing  Equations  45  through  47  with  Equations  58,  61  and  64,  shows 
the  following  equalities: 


c  f  R  .J  IQ  J  =  C  IQ..  ] 
w  1  w  a  L a  ij 


6i'V"V  =  cu>V 1 
5  ■' V1  1  -  “C1 


rVll'CaJrCh«-J  “V'C1 


1  Rf-ll  Cj[  Ch^J  CH^  IQ^  ) 


By  analogy  with  Equations  42  and  58,  the  nonstationary  aerodynamic  lift 
expression,  exclusive  of  the  damping  terms  (velocity  dependent)  to  be 
used  is: 

<V  “  -(X)  c  <k>  [CL„  +  CL(3‘S)'J|  {*J  }]  <65> 

Similarly,  the  moment  is: 

<Hi>  *(-  Fv2)cM/j  <<*>  {V  <66> 
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And  the  hinge  moment  is: 


/l  2\ 

Ha 

~  H/3 

ii 

°H“ 

KjI 

hi  +lh. 

[Q‘j  1 

hi 

(67) 


The  aerodynamic  damping  terms  will  not  include  the  effect  of  spanwise 
coupling  of  strips,  but  will  include  the  effect  of  Mach  number  by  the  use 
of  a  proper  pitch  damping  coefficient,  .  The  aerodynamic  damping 
terms  for  each  strip  are  then:  q 


l9  =o 


(68) 


M. 


ipV2  c2C 
2  H  I  M 


Ipv2  )£  C 
2  H  I V  M 


—  <2  +  0.2W  £-/3 

v  5  V 


ct+  0.2W  fi 
5 


(69) 
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/ 1  2  \  C  _ 

—  pV  | —  c 

,2  H  JV  M 


0.2W  -£-<2+0.6  W  -£/? 

O  v  o  v 


0.2W  <2  +  0. 6W  /? 

5  8 


(70) 


3. 1.2. 1.1  Aerodynamic  Center 

With  additional  reference  to  Statement  "CM  of  Paragraph  3.1.2,  several 
assumptions  will  be  made  concerning  the  location  of  the  aerodynamic 
center.  Theory  does  not  adequately  define  this  parameter,  therefore 
the  analysis  will  be  made  by  varying  the  a.c.  location  within  reasonable 
limits.  The  flutter  speeds  determined  will  correspond  to  the  lowest, 
within  the  defined  a.c.  limits.  The  following  assumptions  will  therefore 
be  made: 

A.  The  a.c.  is  located  at  a  constant  percent  chord  on  all  strips. 
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B.  The  rotation  center,  defined  by  X}  on  Figure  5  is  always  located 
the  same  distance  from  the  midchord  as  the  a.c. ,  therefore: 

W1 

This  assumption  is  exact  for  incompressible  flow  when  Xo  =  1/4, 
Xj"3/4  and  also  at  very  high  Mach  numbers,  when  Xq  =  X^  =  1/2. 

3. 1.2.2  Development  -  Empennage 


Aerodynamic  representation  of  the  empennage  is  used  mainly  to  obtain 
the  aerodynamic  contribution  of  the  empennage  on  the  wing  flutter 
characteristics.  In  addition,  without  tail  effects,  the  airplane  would  be 
unstable  in  the  analog  simulation,  since,  in  effect,  the  simulation  is  a 
model  of  the  airplane  in  free  flight.  The  above  aerodynamic  effects  will 
be  simulated  by  a  lift  force  acting  at  the  aerodynamic  centers  of  the 
horizontal  and  vertical  stabilizers.  The  lift  force  therefore  may  be 
written  as 


L..o.-5'v?CLa«rc<k>  <7: 

where  Cl^Is  a  var*able  and  is  the  effective  tail  angle  of  attack  and 
is  equal  to  (a^,-kQ!w)  where  k  is  the  downwash  angle. 

3.1.3  Mass  Representation 

The  mass  distribution  of  the  aircraft  was  broken  up  into  a  number  of 
lumped  masses  as  is  usual  in  aeroelastic  studies  (lumped  parameter 
system).  Figures  7,  8  and  9  show  the  locations  on  the  aircraft  at 
which  masses  were  lumped.  These  masses  were  derived  from  a  zoning 
of  the  wing,  etc.,  into  small  sections,  and  assuming  that  the  mass  was 
concentrated  at  the  center  of  mass  of  the  section.  Also  considered  are 
lumped  mass  moments  of  inertia  wherever  it  was  felt  neglecting  them 
would  invalidate  a  proper  simulation. 


3 . 2  VIBRATION  (MODAL)  ANALYSIS 


Vibration  modes  were  obtained  from  the  analog  simulation  by  connecting 
a  variable  frequency  voltage  source  which  corresponds  to  a  sinusoidal 
velocity  of  fixed  amplitude  to  a  point  in  the  circuit  where  large  displace¬ 
ments  were  expected.  The  frequency  of  the  voltage  source  was  varied 
until  the  criterion  for  the  presence  of  a  normal  mode  was  satisfied  - 
i.e.,  until  the  input  current  to  the  circuit  (force  applied  to  the  airplane) 
was  a  minimum.  At  each  minimum  of  current,  the  frequency  of  vibra¬ 
tion  and  the  voltages  at  all  nodes  were  recorded.  By  appropriate  scale 
factors,  the  mode  shapes  were  determined. 

3.2.1  Vibration  Modes  (Normal) 

For  determining  the  normal  modes  of  vibration,  the  appropriate  equa  ¬ 
tion  is  Equation  9  with  the  aerodynamic  and  gyroscopic  matrices  set 
equal  to  zero.  Equation  9  therefore  reduces  to  an  expression  governing 
forced  vibration  in  a  vacuum  with  structural  damping: 


(72) 


When  the  impressed  sinusoidal  force  Qj  =  O,  c^. 
and  Equation  72  reduces  to 


O  upon  resonance 


which  is  the  standard  equation  for  normal  modes  with  ft  now  equal  to 
l/w2  for  undamped  vibrations  and  equal  to  (1  +  ig)/ w2  for  damped 
vibrations . 


3.2.2  Vibration  Modes  with  Gyroscopic  Effects 


For  determining  the  aircraft  modes  of  vibration  with  gyroscopic  effects, 
the  appropriate  equation  also  is  Equation  9  with  the  aerodynamic  matrix 
set  equal  to  zero.  Following  the  same  procedure  with  the  impressed 
sinusoidal  force  Q.  as  in  Paragraph  3.2.1,  the  equation  finally  reduces 
to  J 


(74) 


which  is  an  expression  for  vibration  modes  with  gyroscopic  effects  with 
ft  equal  to  for  undamped  vibrations,  and  equal  to  (1  +  ig)/u>2  for 
damped  vibrations. 


3.3 


FLUTTER  ANALYSIS 


The  flutter  analysis  of  the  XV-5A  was  programmed  to  investigate 
several  speed-altitude  combinations  that  would  adequately  represent  the 
proposed  flight  boundary.  Flutter  plots  of  damping  versus  :  oeed  were 
constructed  by  setting  up  the  appropriate  circuits  corresponding  to  a 
particular  speed  and  by  imparting  a  transient,  and  then  measuring  the 
transient  decay  rate  of  free  oscillation.  With  those  data,  plus  the 
appropriate  speed,  a  point  on  the  curve  of  structural  damping  versus 
speed  was  obtained.  Several  speeds  were  selected  in  order  to  define 
adequately  the  curve.  The  recorded  transient  decay  represented  the 
least  damped  component  of  free  vibration.  The  oscillations  were  excited 
by  the  application  of  a  pulse  at  a  point  in  the  aerodynamic  circuit  which 
represented  one  component  or  section.  This  point  was  chosen  to  give 
a  maximum  amplitude  for  the  least  damped  component . 

The  overall  flutter  analysis  included  the  effects  of  varying  the 
following: 

1.  Aileron  mass-balance 

2.  Aileron  spring  restraint 

3.  Distribution  of  wing  bending  material 

4.  Torsional  stiffness  of  the  wing  leading  edge  box 

For  determining  the  flutter  speeds,  the  analogous  equation  is  again 
Equation  9  with  the  gyroscopic  matrix  set  equal  to  zero.  When  the 
impressed  sinusoidal  force  Q.  =  O,  q^  =0  and  the  equation  reduces  to 


Equation  75  shows  that  there  will  be  as  many  roots  as  the  number  of 
physical  coordinates  (q  )  used,  to  described  the  motion  of  the  structure. 
The  root  found  from  the  analog,  as  previously  described,  usually  will 
be  the  lowest  damped  root. 
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4.0  DISCUSSION 


Since  a  plane  of  symmetry  exists  on  the  airplane,  it  is  convenient  and 
logical  to  break  down  the  motion  of  the  airplane  into  symmetric  and 
antisymmetric  motions.  Therefore,  in  the  analog  studies  covered  by 
this  report,  vibration  (modal)  and  flutter  are  treated  as  two  separate 
studies  for  the  symmetric  and  antisymmetric  vibration  and  flutter 
analyses .  The  following  sections  present  the  input  and  the  output, 
(vibration  and  flutter).  The  results  of  the  flutter  analysis  may  be  con¬ 
veniently  divided  into  nine  phases.  Each  phase  has  a  particular  emphasis 
on  some  aspect  of  the  problem.  These  phases  are  outlined  as  follows: 

Phase  I  -  Basic  Structural  Evaluation 

A.  Variable  stiffness  configurations 

1.0  Nominal  Stiffness  Configuration  5A 

2.0  Stiffness  Configuration  5B 

B.  Nominal  mass  condition 
1.0  Mass  Condition  1. 0 

C.  Flight  Condition  D 

1.0  M  =  0. 75;  altitude  =  sea  level 
Phase  II  -  Wing  Front  and  Rear  Spar  Flexibility  Studies 

A.  Variable  stiffness  configurations 
1.0  See  Table  2 

B.  Nominal  mass  condition 
1.0  Mass  Condition  1.0 

C.  Flight  Condition  D 


1.0  M  =  0. 75;  altitude  =  sea  level 


BLANK  PAGE 


Phase  III  -  Aileron  Mass-Balance  Studies 


A.  Optimum  stiffness  configuration 

1.0  Resulting  from  Phase  II  studies 

B.  Variable  mass  conditions 
1 . 0  See  Table  6 

C .  Flight  Condition  D 

1.0  M  =  0. 75;  altitude  =  sea  level 

Phase  IV  -  Re-evaluation  of  Nominal  and  Optimum  Stiffness 
Configurations 

A.  Realistic  stiffness  configuration 
1.0  Resulting  from  Phase  II 

B.  Variable  mass  conditions 
1.0  See  Table  6 

C.  Flight  Condition  D 

1.0  M  =  0.75;  altitude  =  sea  level 
Phase  V  -  Flight  Envelope  Evaluation 

A.  Variable  stiffness  configurations 

1.0  Nominal  Stiffness  Configuration  5A 

2.0  Realistic  stiffness  configuration  resulting  from 

Phase  IV  studies 

B.  Optimum  mass  conditions 

1.0  Resulting  from  Phase  IV  studies 
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C.  Variable  flight  conditions 


1.0  See  Table  3 

Phase  VI  -  Wing  Structural  Nose  Box  Evaluation 

A.  Realistic  stiffness  configuration 

1.0  Resulting  from  Phase  IV  studies 

B.  Optimum  mass  condition 

1.0  Resulting  from  Phase  IV  studies 

C.  Flight  Condition  D 

1.0  M  =  0.75;  altitude  =  sea  level 

Phase  VII  -  Aileron  Equivalent  Spring  Restraint  Studies 

« 

A.  Realistic  stiffness  configuration 

1.  0  Resulting  from  Phase  IV  studies 

B.  Optimum  mass  condition 

1.0  Resulting  from  Phase  IV  studies 

C.  Flight  Condition  D 

1.0  M  =  0.75;  altitude  =  sea  level 
Phase  VIII  -  Aircraft  Simulation  Studies 

A.  Variable  stiffness  configuration 

1.0  Nominal  Stiffness  Configuration  5A 

2.0  Realistic  stiffness  configuration  resulting  from 

Phase  IV  studies 

B.  Optimum  mass  conditions 

1.0  Resulting  from  Phase  IV  studies 
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C .  Variable  flight  conditions 


1.0  See  Table  3 

Phase  IX  -  Aerodynamic  Simulation  Studies 

A.  Nominal  stiffness  configuration  5 A 

B.  Optimum  mass  condition 

1.0  Resulting  from  Phase  IV  studies 

C.  Variable  flight  conditions 
1 . 0  See  Table  3 
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4.1 


IDEALIZATION  OF  AIRCRAFT 


As  mentioned  in  Paragraph  3.  I,  the  aircraft  for  analysis  purposes 
must  be  idealized  to  be  amenable  to  analytical  treatment,  yet  be  of 
sufficient  rigor  to  meet  the  demands  of  engineering  accuracy,  while 
consistent  with  the  capability  of  the  analog  computer.  Figure  G  shows 
the  actual  aircraft  structure  to  be  idealized. 


4.1.1  Structural  Representation 

Since  the  airplane  may  be  thought  of  as  a  number  of  components 
assembled  to  form  the  complete  airplane,  e.g.,  wing,  fuselage,  hori¬ 
zontal  and  vertical  stabilizers  and  control  surfaces,  the  following  sec¬ 
tions  will  deal  with  the  components  as  they  were  treated  during  the 
analog  analysis. 

4. 1.1.1  Wing  -  Aileron 

Subject  to  the  preceding  discussions,  the  actual  wing  structure  shown 
in  Figure  6  was  idealized  into  the  structure  shown  in  Figure  7  for  the 
symmetric  analysis,  and  Figures  8  and  9  for  the  antisymmetric  analy¬ 
sis.  Figures  7,  8  and  9  show  the  discrete  components  which  represented 
the  idealized  model.  Spar  I  on  the  outer  wing  panel  is  a  dummy  spar 
assumed  for  the  structural  idealization,  so  as  to  permit  representation 
of  chordwise  bending.  As  such,  it  was  represented  as  a  beam  without 
webs. 


The  basic  structural  configurations  were  designated  as  Stiffness  Con¬ 
figuration  5A  and  5B.  Configuration  5 A  represented  the  structure  in 
a  4g  maneuver  where  very  little  of  the  skin  is  considered  effective  in 
bending.  Configuration  5B  represented  the  structure  in  a  Ig  flight  con¬ 
dition  where  most  of  the  skin  is  considered  effective  in  bending.  All 
other  stiffness  configurations  studied  were  identical  to  Stiffness  Con¬ 
figuration  5A,  except  for  the  front  and  rear  spar  stiffness  inboard  of 
BL  100.38.  Table  2  lists  the  various  stiffness  configurations  studied 
other  than  the  basic  configuration.  The  listing  is  in  terms  of  a  flexi¬ 
bility  factor  which  was  applied  to  the  basic  configuration  r>A  —  i.e. 

—  =  (— )  ,  x  Flexibility  Factor. 

El  vEl'  nominal 


Figure  10  depicts  the  geometry  of  the  wing  as  used  in  the  structural 
idealization  of  the  aircraft.  The  lengths  shown  represent  average 
values  of  the  mid-plane  lengths  of  Stiffness  Configuration  5.  Due  to 
the  small  sweep  of  the  aileron  and  foi  simplification  of  analysis,  the 
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hinge-line  was  assumed  unswept  as  shown  in  the  aileron  geometry 
sketch,  Figure  11.  The  elastic  restraint,  represented  by  the  overall 
control  circuit,  was  assumed  to  act  at  aileron  station  A55.  The  calcu¬ 
lated  spring  rates  which  were  classified  as  the  nominal  values  were 
K^  =  3.761  x  10^  in-#/rad,  symmetric  and  =  2.614  x  10^  in-#/rad, 
antisymmetric.  These  values  used  in  conjunction  with  the  nominal 
aileron  mass  condition  (Mass  Condition  1.0)  gave  16.0  cps  and  13.4  cps 
as  the  aileron  uncoupled  frequency,  symmetric  and  antisymmetric  re¬ 
spectively.  Figures  12  and  13  show  the  variation  in  bending  stiffness  of 
the  front  and  rear  wing  spars  respectively,  for  the  two  basic  stiffness 
configurations  studied  in  the  analysis.  Figure  14  depicts  the  rib  bend¬ 
ing  stiffnesses  for  the  outer  panel.  As  can  be  seen  with  a  comparison 
with  Figures  12  and  13,  which  are  spar  stiffness  plots,  the  chordwise 
stiffness  varies  from  1.5  to  5.0  times  the  spanwise  stiffness,  indicating 
the  possibility  of  neglecting  chordwise  bending.  Figure  15  depicts  the 
variation  in  bending  stiffness  for  the  center  spar  or  ficticious  stringer 
utilized  in  the  structural  idealization.  The  values  shown  are  for  Stiff¬ 
ness  Configuration  5B  wherein  all  the  skin  is  considered  effective  in 
bending.  Figure  16  presents  the  idealized  wing  configuration  as 
developed  by  the  Structures  Group.  This  idealization  was  followed 
closely  in  the  analysis  with  the  exception  of  the  leading  edge  box  in¬ 
board  of  BL  100.38  (Rib  No.  50)  which,  in  this  analysis,  was  further 
idealized  into  an  equivalent  torque  tube  as  shown  in  Figures  7  and  8. 

This  equivalent  torque  tube  represented  the  same  torsional  stiffness  as 
that  of  the  single  cell  box.  The  variation  in  stiffness  is  shown  in 
Figure  17.  Figure  16  also  presents  the  skin,  spar  and  rib  web  average 
thicknesses  for  use  in  the  torque  box  simulation.  Figure  18  presents 
the  stiffness  values  for  the  aileron  which  was  treated  as  a  beam  flexible 
in  bending  and  torsion  with  the  beam  simply  supported  at  the  hinges. 

The  actual  hinge  locations  shown  in  Figure  18  were  moved  to  that 
shown  in  Figure  11  for  analysis  purposes.  The  difference  in  beam 
length  between  hinges  was  of  the  order  of  1.3  inches,  which  is  negligible. 

4. 1.1.2  Fuselage- Empennage 

Figures  19  through  21  depict  the  geometry  of  the  fuselage  and  of  the  hori¬ 
zontal  and  vertical  stabilizers  as  used  in  the  structural  idealization  of  the 
aircraft.  For  purposes  of  analog  representation,  these  components  of  the 
aircraft  were  treated  as  simple  beams  with  assumed  elastic  axes.  The  fuse 
lage  and  empennage  idealization  was  accomplished  in  a  less  sophisticated 
manner  than  the  wing,  due  to  analog  circuit  limitations.  This  was  felt  to  be 
minor  drawback,  since  the  emphasis  was  to  be  wing  flutter  and  vibration 
characteristics,  rather  than  the  complete  aircraft.  Normally, 


empennage  and  wing  dynamic  analyses  can  be  separate  due  to  the  small 
coupling  between  wing  coordinates  and  those  of  the  empennage.  In  con¬ 
junction  with  the  less  sophisticated  representation  of  the  empennage,  the 
swept  horizontal  stabilizer  was  assumed  to  be  unswept  with  the  elastic 
axis  of  the  horizontal  stabilizer  intersecting  the  vertical  stabilizer 
elastic  axis.  With  this  representation,  the  root  chord  plus  the  closure 
rib  chord  at  BL  70. 16  defined  the  planform  as  shown  in  Figure  20. 
Empennage  control  surfaces,  (elevator  and  rudder),  were  not  considered 
for  the  analysis.  Figures  22  through  24  show  the  stiffness  distributions 
for  the  horizontal  and  vertical  stabilizers  respectively,  as  used  in  the 
analysis.  The  torsional  stiffness  (GJ)  plots  of  both  the  horizontal  and 
vertical  stabilizers  reflect  a  flight  condition  wherein  all  the  skin  was 
considered  effective  in  determining  the  stiffness,  whereas  the  bending 
stiffness  (El)  reflects  a  flight  condition  where  only  the  tension  skin  is 
utilized  in  the  computation.  The  incapability  of  the  torsional  and  bend¬ 
ing  stiffnesses  was  felt  to  have  a  negligible  effect  in  determining  the 
wing  flutter  speeds. 

4.1.2  Aerodynamic  Representation 

As  explained  in  Paragraph  3.1.2,  the  aerodynamic  representation  was 
accomplished  through  the  use  of  steady-state  aerodynamic  influence  co¬ 
efficients  modified  by  an  appropriate  lag  function.  Four  flight  conditions, 
(constant  Mach-altitude)  were  chosen  for  analysis  purposes  that  would 
adequately  represent  the  proposed  flight  boundary,  (Figure  25).  These 
flight  conditions  are  tabulated  in  Table  3  with  accompanying  pertinent 
data  related  to  each  condition.  Also  shown  in  Table  3  is  the  lag  function, 
and  lag  function  coefficients  for  each  flight  condition. 

4. 1.2.1  Wing-Ai  le  ron 

Establishment  of  aerodynamic  strips  at  which  the  aerodynamic  loads 
were  to  be  applied  led  to  the  planform  breakdown  as  shown  in  Figures  26 
and  11  for  the  wing  and  aileron  respectively.  Figures  27  through  33 
present  the  matrices  of  aerodynamic  influence  coefficients  as  explained 
by  Equations  45,  46,  and  47.  These  matrices  are  included  herein  more 
as  a  documentation  effort  than  as  an  aid  in  evaluating  the  overall  wing 
flutter  analysis.  The  wing-alone  pitch  damping  coefficients  as 
explained  by  Equations  69  and  70  are  shown  in  Table  3. 
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4. 1.2.2 


Horizontal  and  Vertical  Stabilizer 


In  line  with  the  less  sophisticated  representation  of  the  empennage  due 

to  analog  circuit  limitations,  the  aerodynamic  strips  chosen  for  the 

horizontal  and  vertical  stabilizers  are  shown  in  Figures  20  and  21 

respectively.  Table  3  also  presents  the  aerodynamic  lift  curve  slope 

coefficient  (Ci  )  used  in  the  aerodynamic  simulation  for  the  tail  load  as 
a 

explained  by  Equation  71. 

4.1.3  Mass  Representation 

Only  one  gross  weight  condition  was  studied  during  the  course  of  the 
analysis.  This  was  the  design  gross  weight  of  approximately  9200 
pounds.  In  line  with  the  usual  lumped  parameter  systems  normal  for 
such  studies,  the  breakdown  for  both  the  wing,  fuselage  and  empennage 
are  shown  in  Tables  4  and  5  for  the  symmetric  and  antisymmetric 
idealization  respectively.  Since  mass-balancing  of  the  aileron  was 
studied  as  a  means  of  evaluating  flutter  fixes,  various  aileron  mass 
configurations  were  used.  These  configurations  are  shown  in  Table  6 
with  the  amount  and  location  of  the  balancing  masses. 


4.2 


VIBRATION  (MODAL)  ANALYSIS 


As  explained  in  Paragraph  3.2,  zero  airspeed  vibration  modes  of  the 
airplane  were  obtained  with  and  without  wing  tan  gyroscopic  effects.  In 
order  to  show  relative  magnitudes  of  deflections  at  various  points  on 
the  airplane,  these  mode  shapes  for  the  nominal  stiffness  configuration 
were  plotted  in  pictorial  form. 

4.2.1  Vibration  Modes  (Nominal) 

Figures  34  through  50  depict  in  pictorial  form  the  normal  mode  shapes 
(without  gyroscopic  forces)  for  Stiffness  Configuration  5 A  and  Mass 
Condition  9,  symmetric  and  antisymmetric.  The  mode  shapes  have 
been  normalized  to  unity  at  the  point  of  maximum  amplitude,  and  are 
plotted  as  such,  or  as  one-half  inch  to  the  inch.  Table  7  lists  for 
comparative  purposes,  the  frequencies,  structural  damping  and  the 
panel  point  at  which  the  airplane  was  excited,  as  explained  in  Paragraph 
3.2.  Damping  values  were  measured  from  the  analog  simulation  as  in 
an  actual  ground  shake  test,  since  the  simulation  represented  a  model 
of  the  airplane. 

During  the  flutter  analysis  of  the  airplane,  structural  changes  in  the 
wing  were  studied  (Phase  II,  III  and  IV  studies)  and  Stiffness  Configura¬ 
tion  5P  represents  one  such  configuration  as  is  explained  in  the  follow¬ 
ing  Paragraph  4.3.  The  normal  modes  of  this  configuration  are  not 
shown  pictorially  in  this  report.  However,  Table  8  lists  for  compara¬ 
tive  purposes.  the  frequencies,  damping,  etc.  of  these  modes. 

4.2.2  Vibration  Modes  With  Gyroscopic  Effects 

Figures  51  through  67  depict  also  in  pictorial  form  the  mode  shapes 
with  gyroscopic  forces  for  Stiffness  Configuration  5A  and  Mass  Condi¬ 
tion  9,  both  symmetric  and  antisymmetric.  Table  9  lists  the  fre¬ 
quencies,  structural  damping,  etc.  As  stated  in  Reference  1,  the 
deflections  shown  for  the  modes  with  gyroscopic  forces  are  the  com¬ 
ponent  of  deflection  in  phase  with  the  panel  point  chosen  as  the  point  of 
excitation.  The  phase  angles  between  recorded  deflections  due  to  gyro¬ 
scopic  forces  were  generally  small  (less  than  10°,  except  for  deflections 
of  very  small  amplitude).  The  in-phase  component  of  deflections  were 
recorded  because  the  phase  angles  were  small,  and  because  deflection 
data  including  phase  angles  cannot  be  represented  pictorially. 


Modal  data  for  Stiffness  Configuration  5P  with  gyroscopic  forces  are 
listed  in  Table  10.  Modal  deflection  data  for  this  latter  stiffness 
configuration  both  with  and  without  gyroscopic  forces  are  presented 
in  Reference  1. 


4.3 


FLUTTER  ANALYSIS 


This  section  is  the  flutter  analysis  portion  of  the  CKA  program,  and  is 
divided  into  a  number  or  phases  which  logically  depict  the  dynamic 
studies  of  the  XV-5A  wing.  The  basic  structural  configurations  herein 
are  Stiffness  Configurations  f>A  and  513,  with  Stiffness  Configuration  5A 
being  termed  nominal.  The  nominal  mass  condition  pertains  to  Mass 
Condition  1.0  (See  Table  (>),  and  all  other  mass  conditions  reflect 
changes  in  aileron  mass-balance  only.  Any  changes  in  wing  structural 
mass  due  to  an  increase  and/or  decrease  in  stiffness  were  neglected  in 
the  interests  of  simplicity  of  analysis.  Such  neglect  was  felt  to  be  justi¬ 
fied  in  view  of  the  broad  nature  of  this  preliminary  analysis.  All  analy¬ 
ses  include  the  appropriate  aileron  degrees  of  freedom  such  as  rotation 
about  the  hinge-line,  deflection  in  the  z  direction,  etc. 

In  general,  most  of  the  results  of  the  flutter  analysis  are  plotted  in  the 
conventional  sense,  i.e.  overall  damping  (g)  versus  Mach  number  (M). 

In  addition,  a  solution  parameter  (V/VQ)  is  plotted  where  VG  corresponds 
to  the  Mach  speed  at  the  particular  altitude  (V0  =  Ma,  where  a  is  the 
speed  of  sound  at  altitude).  The  roots  plotted  represent  different 
branches  of  typical  V-g  plots,  in  that  each  group  of  symbols  represents 
the  response  of  the  aircraft  in  a  different  vibration  mode.  Although 
there  are  an  infinite  number  of  aircraft  modes,  only  those  modes  of 
interest  to  the  flutter  analysis  have  been  plotted.  The  first  subscript 
with  each  point  represents  the  frequency  of  oscillation  in  cycles  per 
second.  The  second,  if  any,  represents  some  variable  as  explained 
on  each  plot. 

4.3.1  Phase  1  -  Basic  Structural  Evaluation 

Phase  I  of  the  flutter  analysis  pertained  to  the  evaluation  of  the  basic 
structural  configurations  with  nominal  mass  conditions.  Since  Flight 
Condition  D  represented  a  5%  flutter  margin  point  (see  Figure  25),  this 
condition  was  chosen  for  the  initial  evaluation  of  the  basic  wing.  During 
this  phase  of  the  studies,  simultaneous  studies  of  a  portion  of  Phase  IX 
were  carried  out.  These  were  the  effects  of  aerodynamic  center  shifts 
on  the  flutter  characteristics  of  the  aircraft. 

4.3. 1. 1  Results 

Figures  68  through  71  present  the  symmetric  and  antisymmetric  V-g 
plots  for  both  basic  structural  configurations,  which  are  Stiffness  Con¬ 
figurations  5A  and  5B  with  nominal  Mass  Condition  1.  0  and  nominal  aile¬ 
ron  spring.  In  this  particular  case,  VQ  is  equal  to  (0.  75)  (a),  and  a  is 
the  speed  of  sound  at  sea  level. 
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4.3.2 


Phase  II  -  Wing  Front  and  Rear  Spar  Flexibility  Studies 


During  this  phase  of  the  analysis,  studies  were  directed  to  evaluating 
the  effects  of  reducing  the  stiffness  (increasing  the  flexibility)  and  there¬ 
by  reducing  weight  on  the  flutter  characteristics  of  the  wing.  Since 
Phase  I  studies  showed  a  flutter  velocity  in  an  antisymmetric  mode,  but 
above  15%  Vl»  it  was  decided  to  determine  what  effects  a  stiffness  re¬ 
duction  on  the  front  and  rear  wing  spars  between  BL  0  and  BL  100.38 
would  have  on  this  flutter  velocity.  The  criterion  used  was  that  flutter 
requirements  would  be  met  at  1. 05  Vl  (M  =  0. 75  at  sea  level)  with 
nominal  aileron  (mass-balance).  Once  this  condition  (marginal  stability 
at  this  point)  is  met,  an  attempt  to  raise  the  flutter  speed  by  aileron 
mass-balancing  would  be  ma^e. 

4.3.2. 1  Results 

Figures  75  through  77  present  the  results  of  the  antisymmetric  analysis 
with  Figures  75  and  76  plotted  as  discussed  in  Section  4.3.  Figure  75 
shows  the  effect  of  increasing  the  flexibility  of  the  rear  spar  from  nomi¬ 
nal  to  220%  of  nominal  while  holding  the  flexibility  of  the  front  spar  con¬ 
stant.  Figure  76  shows  the  effect  of  varying  the  flexibility  of  both  front 
and  rear  spars  an  equal  amount  varying  from  85%  of  the  nominal  values 
to  132.5%  of  the  nominal.  Figure  77  shows  the  damping  (g)  for  various 
flexibility  combinations  for  both  front  and  rear  spars  for  V/VQ  ratios 
of  1. 0  and  1.1.  Only  the  roots  of  the  flutter  mode  (f  «  9. 0-12. 0  cps) 
have  been  plotted,  since  this  mode  showed  a  marked  change  in  overall 
damping  (g)  with  stiffness  changes.  Figures  72  through  74  present 
the  equivalent  results  for  the  symmetric  case  as  was  shown  for  the 
antisymmetric. 


4.3.3  Phase  III  -  Aileron  Mass-Balance  Studies 

As  explained  in  Paragraph  4.3.2,  aileron  mass-balancing  was  studied 
to  determine  its  effect  on  the  flutter  mode  and  consequently  the  flutter 
velocity.  From  Phase  II  studies,  Figures  75  and  76,  it  can  be  seen 
that  Stiffness  Configurations  5N  and  5J  fulfill  the  criterion  stated  in 
Paragraph  4.3.2.  However,  since  Stiffness  Configuration  5N  repre¬ 
sented  a  decrease  of  about  55%  (120%  increase  in  flexibility)  in  rear 
spar  stiffness,  while  Stiffness  Configuration  5J  represented  a  decrease 
of  approximately  25%  in  both  front  and  nar  spar  stiffness  (32.5%  in¬ 
crease  in  flexibility),  it  was  decided  to  use  Stiffness  Configuration  5J 
and  by  means  of  aileron  mass-balance,  raise  the  flutter  speed.  This 


was  done  according  to  a  modified  criterion  which  was  sufficient  mass- 
balance,  properly  located,  so  that  neutral  stability  is  achieved  at 
V/VQ  =  1.1  (approximately  1.15  VL). 

4.3.3. 1  Results 

Figure  79  presents  the  results  of  these  studies.  From  this  figure,  it 
can  be  seen  that  Mass  Conditions  4.0  and  7.  0  both  fulfill  the  stated 
criteria.  Figure  78  presents  similar  results  for  the  symmetric  case, 
but  only  for  these  latter  two  mass  conditions.  Other  mass  conditions 
were  not  studied  during  the  symmetric  analysis,  since  it  was  an  anti¬ 
symmetric  flutter  that  was  evident. 

4.3.4  Phase  IV  -  Re-evaluation  of  Nominal  and  Optimum 
Stiffness  Configurations 

From  Phase  II  and  III  studies,  the  optimum  stiffness  configuration  was 
shown  to  be  5J  [  1.325  (l/EInomina^)  ]  with  either  Mass  Condition  4.  0  or 
7.0.  After  due  consideration,  it  was  decided  that  a  25%  decrease  in 
spar  stiffness  was  more  than  could  be  realistically  tolerated.  Therefore, 
a  closer  examination  was  made  in  order  to  reduce  structural  weight  and 
still  maintain  the  specified  15%  margin  on  flutter  speed.  After  due 
consideration,  it  was  decided  to  reduce  the  stiffness  of  both  spars, 
maintaining  an  overall  damping  (g)  of  approximately  0.02,  the  average 
value  of  structural  damping  in  aircraft  structures  at  a  solution  para¬ 
meter  (V/VQ)  of  1. 1  (approximately  1. 15  at  sea  level).  Once  this 
criterion  was  met,  mass-balancing  of  the  aileron  at  the  tip,  (since  the 
flutter  mode  was  sensitive  to  this  type),  would  be  resorted  to  in  raising 
the  flutter  speed  until  neutral  stability  (g  «  0)  was  reached  at  a  solution 
parameter  of  1.1.  Furthermore,  aileron  tip  balancing  would  be  inves¬ 
tigated  for  the  nominal  Stiffness  Configuration  5A,  with  a  view  toward 
optimizing  the  balance  requirements  for  this  configuration. 

4.3.4. 1  Results 

Figure  80  presents  the  overall  damping  (g)  for  various  flexibility  com¬ 
binations  for  both  front  and  rear  spars  for  a  solution  parameter  (V/V0) 
of  1.1.  The  modes  shown  are  antisymmetric,  since  again  this  was 
shown  to  be  the  sense  of  the  flutter  mode.  From  the  plot,  it  can  be 
seen  that  approximately  a  10%  increase  in  front  spar  flexibility  and 
approximately  a  15%  increase  in  rear  spar  flexibility  would  give  the 
desired  results  as  explained  in  Paragraph  4.3.4.  Figure  82  presents 
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the  actual  results  obtained  from  the  study.  Results  have  been  plotted 
for  both  the  realistic  Stiffness  Configuration  5P  and  the  nominal  Stiff¬ 
ness  Configuration  5A,  both  with  tip  balanced  ailerons,  Mass  Condition 
8.0  and  9.0,  respectively.  Stiffness  Configuration  5P  represents  an 
11%  increase  in  front  spar  flexibility  (15%  decrease  in  stiffness) .  The 
slight  discrepancy  is  due  to  the  manner  in  which  the  overall  damping 
was  determined.  Figure  81  presents  the  equivalent  results  for  the 
symmetric  case,  except  that  no  data  was  taken  with  Mass  Condition  1.0 
for  Stiffness  Configuration  5P. 

4.3.5  Phase  V  -  Flight  Envelope  Evaluation 

Having  established  the  final  stiffness  configurations  with  revised  aile¬ 
ron  mass-balance,  the  flight  envelope  was  evaluated  using  the  aero¬ 
dynamics  available  for  the  other  flight  conditions  (see  Table  3).  Since 
the  aerodynamic  center  is  a  function  of  Mach  Number,  this  parameter 
(a  portion  of  Phase  IX  studies)  was  simultaneously  evaluated  for  the  three 
remaining  flight  conditions. 

4.3.5. 1  Results 

Presented  in  Figures  83  through  86  are  the  V-g  plots,  symmetric  sense 
for  Stiffness  Configurations  5 A  and  5P  for  the  four  flight  conditions 
studied  in  the  analysis.  Flight  Conditions  A,  E  and  B  represent  points 
on  the  flight  envelope  which  encompass  a  15%  flutter  margin  on  Mach 
Number  or  limit  dive  speed  at  sea  level,  5,  000'  and  9,500'  altitude, 
respectively  (see  Figure  25).  Flight  Condition  D  represents  a  5% 
flutter  margin  point  at  sea  level.  Presented  on  each  plot  for  the  latter 
three  flight  conditions  are  the  effects  of  aerodynamic  center  shifts  on 
the  overall  damping  (g).  Figure  87  presents  a  composite  plot  of  the 
four  flight  conditions' V-g  plots  for  an  aerodynamic  center  of  0.20c. 
Figures  88  through  92  present  equivalent  results  for  the  antisymmetric 
analysis.  On  all  plots,  only  the  root  or  mode  of  concern  has  been 
plotted . 

4.3.6  Phase  VI  -  Wing  Structural  Nose  Box  Evaluation 

During  the  course  of  the  analysis,  design  problems  arose  on  a  needed 
cutout  on  the  leading  edge  structural  box  on  the  wing  inner  panel.  Con¬ 
sequently,  it  was  decided  to  investigate  the  effect  of  reduced  stiffness  on 
the  equivalent  torque  tube  (see  Figure  7)  lying  between  BL  61.00  and 
BL  100.38.  Since  a  flutter  point  occurred  during  the  analysis  of  Flight 
Condition  D  (M  =  0.75  at  sea  level),  this  condition  was  chosen  for  study. 


4.3.6. 1  Results 


Figures  93  and  94  present  the  results  of  reducing  the  nose  box  stiffness 
in  the  usual  V-g  plot,  symmetric  and  antisymmetric,  respectively.  The 
factors  shown  are  percentages  of  the  nominal  torsional  stiffness  of  the 
equivalent  torque  tube  given  in  Figure  17.  The  stiffness  configuration 
chosen,  as  basic  for  this  particular  study,  had  no  relationship  with  the 
problem,  other  than  using  one  or  the  other  (i.e.  Stiffness  Configuration 
5A  or  5P,  as  the  base  configuration  upon  which  to  perform  the  study). 

4.3.7  Phase  VII  -  Aileron  Equivalent  Spring  Restraint  Studies 

Since  in  most  cases  of  basic  surface  flutter,  some  degree  of  control  sur¬ 
face  coupling  with  basic  surface  motion  is  involved,  it  was  decided  to 
investigate  the  effects  of  a  variable  aileron  spring  restraint.  As  men¬ 
tioned  in  Paragraph  4. 1,  the  aileron  restraint,  (control  circuit)  was 
idealized  into  an  equivalent  spring  acting  at  the  aileron  hinge  line. 

Again,  as  during  the  Phase  VI  studies,  Flight  Condition  D  (M  =  0.75  at 
sea  level)  wa6  chosen  for  study. 

4.3 . 7. 1  Results 

Plotted  in  Figures  95  and  96  are  the  results  of  varying  the  overall  spring 
rate(K/3),  symmetric  and  antisymmetric,  respectively.  The  factors 
shown  are  various  percentages  of  the  nominal  spring  rates  given  in 
Paragraph  4. 1. 1. 1.  Stiffness  Configuration  5P  was  used  in  conjunction 
with  these  studies  as  in  Phase  VI  studies  providing  a  base  configuration. 

4.3.8  Phase  VIII  -  Aircraft  Simulation  Studies 

Since  a  case  of  wing  flutter  occurred  during  the  course  of  analysis,  it 
was  decided  to  investigate  the  fuselage  degrees  of  freedom  involved. 

Since  a  free-free  airplane  was  simulated  on  the  analog,  it  was  felt 
investigations  along  the  line  of  aircraft  simulation  would  aid  in  providing 
insight  into  the  instability,  and  furthermore  provide  direction  toward 
future  wind  tunnel  flutter  model  simulation.  Studies  covered  the  range 
whereby  the  fuselage  presented  an  infinite  restraint  to  the  wing  (canti¬ 
levered  wing)  to  a  fuselage  and  empennage  free  in  space,  but  having 
infinite  stiffness  in  addition  to  the  nominal  stiffness  present  in  the  fuse¬ 
lage  and  empennage  (see  Figures  22,  23  and  24).  All  flight  conditions 
were  covered  with  Stiffness  Configuration  5P,  providing  the  base  config¬ 
uration  with  the  exception  of  the  cantilevered  case  which  included  Stiff¬ 
ness  Configuration  5A  and  which  was  restricted  to  Flight  Condition  D 
(M  =  0.  75  at  sea  level). 


4.3.8. 1  Results 


Presented  in  Figures  97  through  101  are  the  resulting  V-g  plots  for  the 
aircraft  simulation  studies,  both  symmetric  and  antisymmetric  analyses. 
Figure  97  gives  the  comparison  between  a  cantilevered  wing  and  a  free- 
free  aircraft  with  flexible  fuselage  and  empennage  for  both  Stiffness 
Configurations  5A  and  5P  for  Flight  Condition  D  (M  =  0.75  at  sea  level). 
Figures  98  through  101  presert  the  comparisons  between  a  free-free 
aircraft  with  both  infinite  and  nominal  fuselage  and  empennage  stiffness 
for  all  four  flight  conditions  studied. 

4.3.9  Phase  IX  -  Aerodynamic  Simulation  Studies 

As  stated  in  Paragraph  3.1.2,  simulation  of  the  unsteady  aerodynamics 
required  for  a  flutter  analysis  requires  the  use  of  assumptions  or  empir¬ 
ical  approximations.  The  aerodynamic  center  is  one  such  parameter 
which  is  not  known  adequately  from  theory.  Therefore,  this  parameter 
was  varied  throughout  the  analysis.  Secondly,  the  lag  function  C(k)  (see 
Equation  65  and  Table  1)  is  of  importance  in  analog  simulation  for  flutter, 
and  thus  was  varied  in  the  analysis  from  in  and  out  of  the  circuit. 

Thirdly,  the  pitch  damping  coefficient  (Cj^  ),  wing-alone,  was  not  known 

too  well  for  the  XV-5A  wing  at  the  time.  Studies  were  included  in  the 
analysis  which  also  varied  this  parameter.  A  full  discussion  of  the  pre¬ 
ceding  items  is  presented  in  Reference  7. 

4.3.9. 1  Results 


Variations  in  aerodynamic  center  locations  have  been  shown  throughout 
the  preceding  results  and  will  not  be  duplicated  here.  Figures  102 
through  105  present  the  results  of  removing  the  lag  function  from  the 
aerodynamic  circuits  for  the  four  flight  conditions.  Figure  106  presents 
the  effect  of  varying  the  wing  alone  pitch  damping  coefficient  (CM  )  for 

Flight  Condition  D  (M  =  0.  75  at  sea  level).  A  value  of  tt/8  per  rad.  for 
this  coefficient  is  the  theoretical  value  at  M  =  0.  Figures  107  through 
109  present  similar  results  for  the  remaining  three  flight  conditions, 
with  the  coefficient  being  varied  from  0  to  7r/8,  wing  alone  for  Flight 
Condition  A  (M  =  0.82  at  sea  level);  7r/ 8,  wing  alone  for  Flight  Condition 
E  (M  =  0.90  at  5,000'),  and  finally  zero,  wing  plus  tail  for  Flight  Con¬ 
dition  B  (M  =  0.98  at  9,  500').  The  basic  curves  shown  on  each  plot  of 
Figures  106  through  109  include  the  values  of  the  pitch  damping  coeffi¬ 
cient  (Cj^  ),  wing  alone  listed  in  Table  3. 


5.0  CONCLUSIONS 


Referring  to  Figure  25,  the  proposed  flight  envelope,  it  is  to  be  noted 
that  during  the  course  of  the  program  at  CEA,  a  revised  flight  envelope 
came  into  being.  The  four  flight  conditions  used  in  the  flutter  analysis 
do  not  include  a  portion  of  the  revised  envelope,  M  =0.90  at  12,500' 
altitude  (M  =  1.  035  for  1.15  on  V^)-  However,  it  is  felt  that  the  analysis 
does  adequately  represent  the  revised  envelope,  except  for  the  point  of 
maximum  q  and  Mach  Number,  since  Flight  Condition  B  (M  =  0. 98  at 
9,  500’)  represents  the  same  q  with  the  Mach  Number  difference  being 
0.055. 


5. 1  VIBRATION  (MODAL)  ANALYSIS 

Figures  34  through  67  had  presented  the  modes  for  Stiffness  Configuration 
5A  and  Mass  Condition  1.0,  with  and  without  gyroscopic  forces.  These 
modes  include  the  effect  of  a  flexible  empennage,  which  in  the  CEA 
analysis  was  represented  somewhat  crudely  due  to  equipment  limitations. 
It  is  felt  that  the  modes  which  have  a  large  contribution  of  empennage 
motion  should  be  viewed  somewhat  critically,  and  not  taken  as  being 
representative  of  the  actual  structure.  Since  this  analysis  is  concerned 
with  wing  characteristics,  for  which  analog  representation  was  quite 
thorough,  emphasis  must  be  placed  on  wing  results. 

5.1.1  Vibration  Modes  (Normal) 

The  normal  modes  of  concern  in  the  analysis  are  those  shown  in  Figures 
35  and  38  for  the  symmetric  case,  and  those  shown  in  Figures  44  and  46 
for  the  antisymmetric  case,  all  for  Stiffness  Configuration  5A  and  Mass 
Condition  9.0.  Figure  35,  frequency  8.50  cycles  per  second,  depicts 
the  mode  which  might  be  termed  first  wing  bending.  Figure  38,  fre¬ 
quency  16.55  cycles  per  second,  depicts  the  mode  which  might  be  termed 
first  wing  torsion.  Similarly,  for  the  antisymmetric  ca6e,  Figure  44, 
frequency  11.66  cycles  per  second;  and  Figure  46,  frequency  16.55 
cycles  per  second,  depicts  what  might  be  termed  the  first  wing  bending 
and  torsional  modes,  respectively.  Some  of  the  other  modes  shown  are 
also  wing  modes,  but  of  a  higher  frequency  with  the  remaining  exhibiting 
large  amounts  of  coupling  with  the  other  components,  i.e.  aileron,  fuse¬ 
lage  and  empennage.  Aileron  modes  which  should  be  mentioned  are  those 
of  Figure  37,  frequency  12.38  cycles  per  second,  and  Figure  43,  fre¬ 
quency  10.34  cycles  per  second,  symmetric  and  antisymmetric,  res¬ 
pectively.  The  above  modes  have  been  emphasized,  since  in  normal 


cases,  it  is  these  modes  that  comprise  the  mode  of  flutter,  as  will  be 
shown  in  subsequent  sections. 

5.1.2  Vibration  Modes  with  Gyroscopic  Effects 

Figure  52,  frequency  8.42  cycles  per  second  and  Figure  55,  frequency 
16.34  cycles  per  second,  are  the  important  modes  of  concern  for  the 
symmetric  case,  whereas  Figure  61,  frequency  11.65  cycles  per  second 
and  Figure  63,  frequency  16.32  cycles  per  second  depict  the  equivalent 
antisymmetric  wing  modes,  again  for  Stiffness  Configuration  5 A  and  Mass 
Condition  9.  0.  The  aileron  modes  are  shown  in  Figure  54,  frequency 
12.40  cycles  per  second  and  Figure  60,  frequency  10.32  cycles  per 
second,  symmetric  and  antisymmetric,  respectively.  From  a  compar¬ 
ison  of  the  above  modes  with  the  modes  discussed  in  Paragraph  5.1.1, 
it  is  apparent  that  gyroscopic  forces  are  not  too  significant  in  causing 
radical  modal  changes.  This  is  so,  since  the  gyroscopic  moments  are 
small.  The  effective  mass  moment  of  inertia  of  the  fan  rotating  parts 
about  the  fan  rotor  axis  is  30.20  slug-ft^  and  the  fan  speed  for  100% 
power  is  2640  rpm. 

5.2  FLUTTER  ANALYSIS 


To  adequately  evaluate  the  flutter  characteristics  of  the  XV-5A  aircraft 
from  the  overall  CEA  analysis,  it  is  best  to  follow  the  results  of  each 
phase  of  the  flutter  analysis  as  presented  in  Paragraph  4.3. 

5.2.1  Phase  I  -  Basic  Structural  Evaluation 


Figures  68  through  71  depict  the  V-g  (or  M-g)  plots  for  the  basic  stiff¬ 
ness  configurations  of  the  XV-5A  aircraft  for  various  aerodynamic  cen¬ 
ter  (a.c.)  locations  for  Flight  Condition  D  (M  =  0.75  at  sea  level).  The 
range  of  a.c.  locations,  0.20c  to  0.30c  differ  slightly  from  that  pre¬ 
sented  in  Figure  30  (CM  I  ]  ),  in  that  it  was  felt  the  wider  range 


in  a.c.  travel  would  better  define  the  dependence,  if  any,  of  the  flutter 
speed  on  the  a.c.  travel.  However,  as  Figures  68  and  69  show,  a.c. 
location  is  insignificant  in  altering  the  trends  of  the  symmetric  case. 

This  is  not  the  case  in  the  antisymmetric  analysis  as  shown  in  Figures 
70  and  71,  in  which  a  flutter  point  does  occur,  indicating  that  a.c.  travel 
is  significant  when  a  forward  location,  0.20c  produces  the  lowest  flutter 
speed  (extrapolated  in  Figure  71).  The  difference  in  speeds  between  the 
most  forward  and  the  most  aft  location  is  small,  yet  this  study  does  point 
out  the  need  for  correctly  estimating  the  aerodynamic  center  location. 

In  addition,  as  Figure  70  shows,  the  more  critical  stiffness  configuration 


Is  5A  (4g  representation).  It  is  this  configuration  that  is  used  as  the 
nominal  configuration  throughout  the  entire  analysis.  The  flutter  mode 
is  comprised  of  antisymmetric  wing  bending  coupled  with  aileron  rota¬ 
tion.  Although  the  flutter  frequency  is  about  11.6  cycles  per  second,  and 
the  antisymmetric  wing  mode  was  shown  to  be  11.66  cycles  per  second 
with  antisymmetric  aileron  10.34  cycles  per  second,  it  must  be  remem¬ 
bered  that  the  predominant  motion  is  as  mentioned,  which  implies  coupling 
with  aileron,  fuselage  and  empennage.  The  flutter  speed  shown  in  Fig¬ 
ure  70,  M  =  0.  845  does  not  represent  the  true  flutter  speed,  since  M  = 
0.75  aerodynamics  were  used  in  the  analysis.  Therefore,  additional 
Mach  Numbers  must  be  run.  The  results  when  plotted  on  an  auxiliary 
plot,  as  in  Figure  110,  will  define  the  true  flutter  speed,  if  any,  at  sea 
level.  Conclusions  reached  on  a  true  flutter  speed  are  discussed  in 
Paragraph  5.2.5. 

5.2.2  Phase  II  -  Wing  Front  and  Rear  Spar  Flexibilities  Studies 

Phase  I  studies  have  shown  that  a  flutter  margin  in  excess  of  15%  on 
does  exist  and  therefore  a  weight  reduction  study  was  instigated.  Phase 
I  had  shown  that  the  critical  mode  of  flutter  was  an  antisymmetric  bend¬ 
ing  coupled  with  aileron  rotation.  Obviously,  a  reduction  in  bending 
stiffness  would  lower  the  flutter  speed.  As  mentioned  in  Paragraph 
4.3.2,  the  criterion  was  to  have  neutral  stability  at  1.05  VL  (M  =  0.  75 
at  sea  level),  this  now  being  a  true  flutter  speed,  and  then  attempt  to 
raise  the  flutter  speed  by  proper  aileron  mass-balancing.  The  signifi¬ 
cant  reduction  in  flutter  speed  is  shown  in  Figures  75  and  76.  The 
flutter  frequency  reduced  to  9.5  cycles  per  second  for  a  120%  reduction 
in  rear  spar  flexibility,  and  10.5  cycles  per  second  for  a  32.5%  reduc¬ 
tion  in  both  front  and  rear  spar  flexibilities.  The  symmetric  results 
are  not  altered  appreciably  as  shown  by  Figures  72  through  74. 

5.2.3  Phase  III  -  Aileron  Mass- Balance  Studies 

Figure  79  shows  the  results  of  applying  proper  aileron  mass-balance  to 
Stiffness  Configuration  5J  (32.5%  increase  in  front  and  rear  spar  flexi¬ 
bilities).  These  results  show  the  definite  inertial  coupling  between  the 
wing  bending  mode  and  the  aileron  rotational  mode.  Mass  Condition 
4.0  (see  Table  6.0)  fulfills  the  criterion  stated  in  Paragraph  4.3.3,  and 
also  points  out  the  established  fact  that  for  a  flutter  mode  involving 
coupling  between  main  surface  bending  and  control  surface  rotation, 
mass-balance  at  the  outboard  tip  of  the  control  is  the  optimum  position. 
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5.2.4  Phase  IV  -  Re-evaluation  of  Nominal  and  Optimum  Stiffness 
Configurations 

Phase  II  and  III  studies  have  shown  that  a  stiffness  reduction  and  conse¬ 
quently  a  weight  savings  was  feasible.  These  studies  are  mainly  aca¬ 
demic  in  that  a  32.5%  increase  in  spar  flexibilities  (1/EI)  or  a  25% 
reduction  in  I  (E  being  constant)  is  more  than  could  be  tolerated  from  a 
practicality  standpoint.  However,  these  studies  do  indicate  that  some 
reduction  in  stiffness  is  feasible.  Figure  80  is  a  plot  of  damping  versus 
front  spar  flexibility  factor,  and  it  shows  that  a  small  increase  in  flex¬ 
ibility  within  the  range  of  practicality,  is  feasible  for  neutral  stability  with 
adequate  aileron  mass-balance  at  1. 15  (1.15  M^).  Stiffness  Config¬ 

uration  5P  with  Mass  Condition  8.0,  as  shown  in  Figure  82,  fulfills  the 
criterion  stated  in  Paragraph  4.3.4.  The  11%  front  spar  and  18%  rear 
spar  flexibility  increase  lepresents,  respectively,  a  10%  and  15% 
decrease  in  I  for  the  front  and  rear  spars  (inboard  panel).  The  weight 
savings  thus  gained  represent  the  weight  of  the  spar  reduction  (in  moment 
of  inertia),  plus  the  difference  in  total  mass-balance  per  aileron  (7.  0  lb. 
for  Mass  Condition  8.0  versus  11.08  lb.  for  Mass  Condition  1.0).  Fig¬ 
ure  82  also  depicts  the  results  of  applying  the  information  gained  during 
the  Phase  III  studies  to  the  nominal  stiffness  configuration  (5A).  Here, 
proper  mass-balancing  of  the  aileron  resulted  in  a  weight  savings  of 
6.08  pounds  per  aileron,  with  a  slight  reduction  in  flutter  speed,  but 
still  at  15%  above  as  shown  in  Figure  82. 

5.2.5  Phase  V  -  Flight  Envelope  Evaluation 

Having  established  the  optimum  stiffness  configuration  and  mass  condi¬ 
tion,  an  analysis  was  performed  for  all  flight  conditions.  Figures  87 
and  92  show  the  results  of  the  symmetric  and  antisymmetric  analyses, 
respectively,  for  both  the  optimum  Stiffness  Configuration  5P  and  the 
nominal  Stiffness  Configuration  5A,  both  with  aileron  tip  mass-balance. 
Figure  92  depicts  the  trend  of  the  antisymmetric  flutter  mode.  A  flutter 
point  occurs  only  at  Flight  Condition  D  (M  =  0.75  at  sea  level). 

All  other  flight  conditions  are  free  of  flutter,  with  no  flutter  occurring 
during  the  symmetric  analysis  for  all  flight  conditions  as  depicted  in  Fig¬ 
ure  87.  Plotting  of  the  flutter  point,  occurring  during  Flight  Condition  D 
(M  =  0.75  at  sea  level),  on  an  auxiliary  plot  as  in  Figure  110,  shows  that 
no  true  flutter  speed  exists  for  either  stiffness  configuration  with  aileron 
tip  mass-balance.  A  true  flutter  speed  would  occur,  for  example,  if 
neutral  stability  was  reached  at  a  flight  parameter  of  V/VG  =  1.0,  for 
example,  when  M  =0.75  at  sea  level.  Figure  79  shows  a  flutter  point 


occurring  at  V/VQ  =  1.0  for  M  =  0. 75  at  sea  level.  This  was  shown  to  be 
a  flutter  point  corresponding  to  a  stiffness  configuration  with  increased 
spar  flexibilities  (Stiffness  Configuration  5J)  and  Mass  Condition  1.0. 

This  point  is  plotted  in  Figure  110  to  demonstrate  the  meaning  of  this 
auxiliary  plot. 

5.2.6  Phase  VI  -  Wing  Structural  Nose  Box  Evaluation 

Studies  of  the  reduction  in  the  wing  structural  nose  box,  Figure  94,  have 
shown  that  a  slight  decrease  in  flutter  speed  occurs  in  going  from  nomi¬ 
nal  torsional  stiffness  to  zero  torsional  stiffness.  This  is  so,  since  the 
flutter  mode  is  one  of  wing  bending  coupled  with  aileron  rotation  and  the 
nose  box  contributes  more  to  the  overall  torsional  stiffness  of  the  wing. 
The  torsional  stiffness  of  the  wing  inner  panel  is  generated  mainly  by 
differential  bending  of  both  front  and  rear  spars,  and  reductions  in  spar 
stiffness  would  affect  the  torsional  mode  (see  Tables  7  and  8  for  modal 
frequency  comparisons).  However,  Phase  II  studies  (see  Figures  75  and 
76)  have  shown  that  the  wing  torsional  mode  is  not  critical,  even  with 
such  reductions. 

5.2.7  Phase  VII  -  Aileron  Equivalent  Spring  Restraint  Studies 

Phase  VII  investigated  the  effects  of  varying  the  aileron  spring  restraint, 
since  in  essence  this  value  is  somewhat  difficult  to  calculate  accurately. 
Figure  96  depicts  the  variation  in  flutter  speed  obtained  by  reducing,  and 
also  increasing,  the  overall  spring  constant  (K^).  Also  shown  are  the 
effects  of  having  a  free-floating  aileron  and  a  rigidly  fixed  aileron.  The 
fhitter  frequencies  indicate  a  cha.ige  in  flutter  mode  from  the  basic  wing 
bending-aileron  rotation  mode,  to  one  involving  wing  bending  coupled 
with  another  mode  in  going  from  a  free-floating  aileron  to  a  rigidly 
fixed  aileron.  The  trend  shown  in  Figure  96,  (that  of  a  minimum  flutter 
speed  for  some  factor  of  the  nominal  spring  constant  (K^)  and  then  in¬ 
creasing  flutter  speeds  above  and  below  this  minimum),  is  common  when 
dealing  with  basic  surface-control  surface  flutter. 

5.2.8  Phase  VIII  -  Aircraft  Simulation  Studies 


This  phase  of  the  overall  analysis  indicates  that  a  proper  aircraft  simu¬ 
lation  must  be  made  in  order  to  determine  the  contribution  of  component 
degrees  of  freedom  that  go  into  making  up  the  mode  of  flutter.  Figure  97 
adequately  demonstrates  this.  By  simply  cantilevering  the  wing,  flutter 
does  not  exist,  indicating  that  fuselage  or  aircraft  degrees  of  freedom 
contribute  to  the  flutter  mode.  Furthermore,  Figure  98  indicates  that  by 
simply  including  aircraft  rigid  body  degrees  of  freedom  (side  translation, 
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yaw  and  roll  for  the  antisymmetric  case),  the  mode  of  flutter  is  preserved, 
with  a  slight  decrease  in  flutter  speed  due  to  rigidizing  of  the  fuselage 
and  empennage. 

5.2.9  Phase  DC  -  Aerodynamic  Simulation  Studies 

This  section  is  presented  more  as  a  discussion  on  the  importance  of 
properly  simulating  certain  terms  in  the  aerodynamic  representation, 
rather  than  as  a  discussion  on  the  flutter  characteristics  of  the  XV-5A 
aircraft. 

Effects  of  aerodynamic  center  shifts  have  been  discussed  in  Paragraph 
5.2. 1,  with  the  more  critical  a.c.  location  being  taken  as  0.20c,  as 
shown  in  Figures  70  and  71  for  Flight  Condition  D  (M  =  0. 75  at  sea  level). 
This  location  corresponds  to  the  lowest  flutter  speed  and  is  felt  to  be 
reasonable  for  this  flight  condition.  The  percentage  change  from  that 
obtained  with  the  a.c.  location  at  0.25c,  the  theoretical  subsonic 

value,  is  very  small,  yet  this  study  does  indicate  a  need  for  correct 
simulation  of  this  parameter.  Figures  84  through  86  present  the  effects 
of  a.c.  location  in  the  symmetric  sense  for  the  remaining  three  flight 
conditions,  whereas  similar  results  for  the  antisymmetric  analyses  are 
presented  in  Figures  89  through  91.  Here,  the  aerodynamic  center 
locations  range  from  0.20c  to  0.32c.  The  aft  location  is  more  indica¬ 
tive  of  the  trend  with  increasing  Mach  Number. 

Figure  102  demonstrates  the  importance  of  adequately  representing  the 
lag  function  when  undergoing  an  analog  flutter  analysis.  With  the  lag 
function  completely  removed  from  the  circuit,  the  flutter  speed  is 
drastically  reduced  in  the  antisymmetric  case,  and  a  flutter  point  occurs 
(extrapolated)  in  the  symmetric  case  for  Flight  Condition  (M  =  0.  75  at 
sea  level).  The  symmetric  case  appears  to  be  in  error  when  consider¬ 
ing  the  change  in  flutter  frequency  from  results  obtained  with  and  without 
the  lag  function,  when  compared  with  the  results  obtained  for  the  other 
flight  conditions,  Figures  103  through  105.  The  overall  trend,  though, 
is  to  diminish  the  effect  of  the  lag  function  as  the  Mach  Number  is  in¬ 
creased  as  mentioned  in  Reference  1.0. 

The  third  aerodynamic  parameter  to  be  varied,  that  of  the  pitch  damping 
coefficient  (CM  ),  shows  the  effect  of  reducing  this  value  from  that  used 

q 

in  the  analysis  (see  Table  3)  to  tt/8,  the  theoretical  value  at  M  =  0,  and 
finally  reducing  it  to  zero.  As  reference  7,  Figure  11,  depicts,  the 
value  of  CM  increases  from  tt/8  at  M  =  0  to  infinity  as  M  =  1.0  is 

q 


approached.  The  simulation  of  a  proper  pitch  damping  coefficient  is 
apparent  when  one  examines  Figures  106  through  109  in  which  flutter 
results  are  plotted  for  four  flight  conditions.  Taking  into  account  the 
trend  of  C™  as  M  =  1. 0  is  approached,  and  the  existence  of  a  flutter 

q 

trend  as  shown  in  Figures  106  and  107,  it  can  be  seen  that  correct  simu¬ 
lation  of  this  coefficient  is  apparent.  With  the  reduction  of  the  value  of 
the  coefficient  to  nr/ 16  and  to  zero,  a  second  flutter  mode  emerges  as 
depicted  in  Figure  107,  antisymmetric  with  the  critical  mode  showing  a 
marked  decrease  in  flutter  speed.  In  addition,  from  Figure  107  with  the 
inclusion  of  reduced  values  of  C-.  ,  a  flutter  mode  appears  for  a  con- 

q 

figuration  which  is  essentially  stable  in  the  symmetric  sense. 

5.2.10  Overall  Evaluation 

The  results  of  Phase  I  through  IX  have  indicated  that  the  XV-5A  aircraft 
(wing),  on  the  basis  of  the  preliminary  data  utilized  in  the  analysis  is 
free  of  flutter  within  the  proposed  flight  envelope  and  revised  envelope 
as  shown  in  Figure  110.  Studies  have  shown  that  mass-balancing  of  the 
aileron  is  a  critical  factor  and  must  be  carefully  evaluated  for  future 
design  efforts  on  the  XV-5A  aircraft.  Similarly,  the  overall  spring 
restraint  (K^)  does  affect  the  flutter  characteristics,  but  the  variation, 
at  most,  will  result  in  flutter  speeds  above  15%  on  VL  of  ML.  The 
flutter  mode,  that  of  antisymmetric  wing  bending,  coupled  with  aileron 
rotation,  indicates  that  proper  aircraft  simulation  must  be  used  in  future 
flutter  evaluation  of  the  XV-5A  aircraft  with  due  regard  given  to  fuse¬ 
lage  and/or  aircraft  degrees  of  freedom. 
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Figure  2  Beam  Element 


Figure  3  Torque  Tube  Element 
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Figure  13  Rear  Wing  Spar  -  Bending  Stiffness  Distribution 
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Figure  34  Symmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces 
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Figure  39  Symmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  Out  -  Mode  6 
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Figure  42  Antisymmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  Out  -  Mode  1 


Figure  44  Antisymmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  Out  -  Mode  3 


Figure  46  Antisymmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  Out  -  Mode  5 
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Figure  47  Antisymmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  Out  -  Mode  6 
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Figure  48  Antisymmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  Out  -  Mode  7 
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Figure  52  Symmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  In  -  Mode  2 


Figure  53  Symmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  In  -  Mode  3 
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Figure  54  Symmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  In  -  Mode  4 
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Figure  55  Symmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  In  -  Mode  5 


Figure  56  Symmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  In  -  Mode  6 


Figure  57  Symmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  In  -  Mode  7 
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Hill 


Figure  62  Antisymmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  In  -  Mode  4 
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Figure  63  Antisymmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  In  -  Mode  5 
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Figure  65  Antisymmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  In 


Mill 


Mode  No 


Figure  67  Antisymmetric  Mode  Shapes  -  Stiffness  Configuration  5A  -  Gyroscopic  Forces  In  -  Mode  9 
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Figure  70  Phase  I  Studies  -  Antisymmetric:  Flight  Condition  D  (M  =  0.75  at  Sea  Level) 


Figure  71  Phase  I  Studies  -  Antisymmetric 


Flight  Condition  D  (M  =  0.75  at  Sea  Level) 
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Figure  74  Phase  II  Studies  -  Symmetric:  Flight  Condition  D  (M  =  0.75  at  Sea  Level) 
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Figure  77  Phase  II  Studies  -  Antisymmetric:  Flight  Condition  D  (M  =  0.75  at  Sea  Level) 
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Figure  81  Phase  IV  Studies  -  Symmetric:  Flight  Condition  D  (M  =  0.75  at  Sea  Level) 
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Figure  84  Phase  V  Studies  -  Symmetric:  Flight  Condition  A  (M  =  0.82  at  Sea  Level) 


Figure  85  Phase  V  Studies  -  Symmetric:  Flight  Condition  E  (M  =  0.90  at  5,000' 
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night  Condition  B  (M  =  0. 98  at  9, 500') 


Figure  90  Phase  V  Studies  -  Antisymmetric:  Flight  Condition  E  (M  =  0.90  at  5,000') 


Figure  95  Phase  VII  Studies  -  Symmetric;  Flight  Condition  D  (M  =  0.75  at  Sea  Level) 
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Figure  107  Phase  IX  Studies  -  C 


-  Flight  Condition  A  (M  =  0. 82  at  Sea  Level) 


TABLE  1 


AERODYNAMIC  DATA 


X 

a 

W2/Wl 

0. 8  W  ,W 

5/  4 

0.4W8/W7 

W6 

0.  2  Wc 

5 

0.65 

0.241 

0.252 

0.250 

0.  13 

0.206 

0.75 

0. 173 

0.163 

0. 160 

0. 137 

0. 122 

0.  116 

TABLE  2 

STIFFNESS  CONFIGURATIONS 

Configuration  Number  Flexibility  Factor 

Front  Spar  Rear  Spar 


A 

Basic  (Nominal) 

Basic 

bB 

Basic 

Ban  i.e 

bC 

0.8b 

i-ib 

A 

l-ib 

1.1b 

l.Lb 

o.8b 

bF 

o.8b 

0.8b 

iXJ 

1.00 

l.Lb 

1.40 

1 .  ho 

'A 

i-3b 

1.3b 

A 

i.32b 

i.32b 

bK 

1.00 

l.so 

% 

1.00 

1.90 

[M 

1.00 

h .  10 

bN 

1.00 

L» .  20 

bP 

1.11 

1.1.8 

1 


x  Flex.  Factor 


El 


JL 

El 


Nomina L 


TABLE  3 


FLIGHT  CONDITIO] 


Flight 

Condition 

Mach 

No. 

Altitude 

(Ft.) 

Vo 

(Knots) 

( 

% 

(lb/Ft2) 

D 

0.75 

0 

496 

833 

A 

0.82 

0 

542 

996 

E 

584 

996 

B 

0.98 

626 

996 

Lag  Function 
A(^)P  ♦  1 
B(^  )P  +  1 


dt 


'VTCLa<Tall) 

(l/rad.i  (l/racl.) 


Lag  Function 
Coefficient 


A 


B 


TABLE  4 


MASS  AND  INERTIA  DISTRIBUTION 
(SYMMETRIC  IDEALIZATION) 


Wing : 

(Weight:;  Dt 

titr  Ibution 

5/0/O2)  * 

baaie 

Component  ... 

Winy 

Wing 

Fan 

Pune  luge 

A  Uoted  Total  AFan 

Item:;  Wing  (5/7A&) 

Fuue  luge 
Ai  l.o  ted 

I  teirui 

Rev  i  tied 
Total 
Wing 

Pune  1 
Point 

Wt.-Lb. 

Wt.-Lb. 

Wt.-Lb.  Wt.-Lb.  Wt.-Lb. 

Wt.-Lb. 

Wt . -Lb . 

20G 

51.8 

51.8 

201 

0 

0 

2  OK 

19-8 

19.6 

50G 

Bi.O 

81.0 

501 

0 

720.0 

59-9 

799-9 

50K 

74 . 4 

74.4 

1*00 

17-7 

17.7 

40K 

20. 1 

20.1 

50G 

24 . 1 

24.1 

501 

L9.8 

19.  B 

5  OK 

'*7.5 

'*7-5 

60G 

15-9 

15.9 

ooi 

10.9 

16.9 

OOK 

11.1. 

11.1 

YOG 

9.0 

9.6 

701 

11.9 

11.9 

YOK 

9-5 

9-5 

Bog 

0.4 

6.4 

001 

9-8 

5.8 

bok 

6.4 

0.4 

90G 

4 . 9 

4.5 

901 

4.5 

4.5 

90K 

5.9 

5.9 

>: 

459.4 

720.0 

97.5  1290.7  59-9 

-97-5 

1195.5 

Wing 

Fan:  (Lumped  at  Panel  Point  -  See  Above) 

W  =  ’ 

799.9  Lb. 

I  =  0.  9  X  L(P  Lb -In' 

y  ,  (  o 

Ix  =  0.906 5  X  10°  Lb-In^ 
*  Per  Side 


163 


TABLE  4  (Cont.) 

MASS  AND  INERTIA  DISTRIBUTION 
(SYMMETRIC  IDEALIZATION) 


Fuselage 

:  (Weights 

Distribution  4/13/69)  * 

Component 

Basic 

Fuselage 

Wing  Fun 

Flap 

Vertical 
Wing  Stab. 

Total  Fu  se  i  age 

Pane;  l 
Point 

Wt.-Lb. 

Wt.-Lb. 

Wt.-Lb. 

Wt.-Lb.  Wt.-Lb. 

Wt.-Lb. 

iOA 

9%  00 

99.06 

10C 

602.  10 

602. 16 

10E 

9i'*.02 

914.62 

10G 

74b. 64 

47.0 

799.04 

101 

)04.  00 

90 .9 

b.8 

40). yb 

10K 

**79.90 

*90 

41.9 

9}6.)0 

LOM 

90b. 99 

20b. 99 

100 

li‘9.07 

129.07 

10Q 

)4.00 

9b. 2 

92.20 

£ 

910b. 48 

90. ) 

iy.} 

970  98.2 

Horizontal. 

Stabilizer: 

(Weight:; 

Distribution  4/17/62)  * 

Panel  Point 

Wt.-Lb. 

I  -Lb. -In 

ea 

H£  .10 

-jy.y 

by;'.’ 

H.S.pQ 

:m  .5 

6024 

L 

si 

I4y4o.o 

Total  Airplane  I)  i  :;t.ri  l-u  t  i  on :  (l/,1  Airplane) 


Fuselage;  5108.6  Lb 

Wing  Kan  to  Fuselage;  90 .  5 

FLnp  to  Fuselage;  19*3 

Wing  Items:  97*3 

Vert  I eul  Stabilizer:  SH.;* 

Total  Fuse  lug.;;  L  5575.6 

Horizon  tat  Stub!  1  i  zer :  61 

All  .Ton  (NomLnul);  i'4.6 


Total  Gross  Weight;  p  (537 5.0  ♦  M.P  4 
Center  *<1‘  GmvLt.y:  F.S.  959.04 


Wing: 

Wing  Fan: 

Wing  I terns ; 
Total  Wing: 

A  Fan ; 

Wing  Items; 
Rev. Total  Wing 


439.».  Lb. 

790.0 

97.3 

L  lL'^i.7 
59.9 
-97  •  3 
i:  1199.3 


l‘4.9  +  1199.5)  =  p  (4644.6)  =  9989.9 


*  Fun.;  1  ug» *  -  1  /;*  Airplane}  Horizontal  stabilizer  -  Per  Side 
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TABLE  5 

MASS  AND  INERTIA  DISTRIBUTION 
(ANTISYMMETRIC  IDEALIZATION) 


Win*':  (W'-ii'lit.:;  liiiilrilmliuii  )  * 


Cwmpunnnt 

Hun  i  c 

W  i  f  m  r. 

Punt:  1  Uf 'o 

A  1  1  ot.i.-ii 

Win#'  Fun  Ilfim; 

Total 

W i  >W'  Win*',  Fan 

Flap 

AFun 

( 9/Y/oo 

Ki-v  i  Lii.'d 
Total 
)  Win« 

Pan*.*  1 
Pol  nt 

Wt .  -Lti . 

Wt.-LH.  Wt.-LO. 

WL.-l.h.  Wl.-U*. 

Wt .  -LI* . 

Wt.-L,H. 

Wt.-J.b. 

KXj 

0 

4‘/\0 

i*y.o 

101 

u 

H.H 

'10.  1 

99 .  i 

IUK 

0 

41.0 

19.1 

00.6 

MUG 

'll  ,M 

11.6 

.'01 

0 

0 

.'OK 

10. t 

19.'' 

'iUj 

Ml  .0 

bl.o 

'>0L 

u 

f.-u.d 

19-9 

799-9 

'>0K 

(U.h 

[h.h 

l  Y  •  Y 

17-7 

4uK 

I 

00.  1 

*'UJ 

' 

i-ii.l 

‘>01 

i  » 

19-6 

OOK 

**Y  * 

'*7.1 

OUi 

1  0.0 

1  1.9 

OUl 

10.0 

10.9 

OoK 

11.1 

11.1 

YOJ 

0.0 

9-0 

Yol 

11.0 

11.9 

YOK 

0. 

9.' 

Mix; 

ojl 

0.1* 

moi 

o  .m 

9.6 

MuK 

(i .  >i 

o.i* 

'MJ 

)|.0 

'*.9 

OOI 

4.5 

'*.1 

V  'K 

5.0 

1.9 

>* 

n^o  ji 

|Vo.o  OY  •  '> 

I.  ■‘A>.  Y  90.1 

19.  .5 

y.j.y 

i402.lj 

W i  1 14 •;  Fun:  ( -<i  at  Puncl  Point 

0 

> 

0 

< 

1 

O 

s-N 

7V;. V  1.0. 
oJt;-iy  x  iu() 


w 
K 

r  0.900 >  X  10  Lh.-In 

X 

*  P«  r  Ji Id*; 


TABLE  5  (Cont.) 

MASS  AND  INERTIA  DISTRIBUTION 
(ANTISYMMETRIC  IDEALIZATION) 


FuucLugo:  (Weighty  Dkit.-ltmtl  on  4/L5/O2)  * 


Component 

Basic 

< 

+ 

< 

Total  Component  Basic  Pune 

Pane  1 

Pune  lu£t 

Pune  La^'e 

Panel 

-b 

Point 

Wt.-Lb. 

Wt.-Lb.  Wt.-Lb. 

Wt.-Lb. 

Point 

1  x  10  -L 

X 

0A 

99-06 

90.00 

1.49.06 

OA 

0.0810 

OC 

602.16 

-  90.00 

992. L6 

OC 

0.4949 

0E 

91.4.62 

914.62 

OE 

0.4899 

0G 

749.64 

L49.07  '.>89.80 

1488.1.1 

OG 

0.9274 

01 

504.88 

-504.88 

0 

01 

0.5450 

OK 

479.90 

199.21  702.46 

1555- 17 

OK 

O.3279 

0M 

208.99 

208.99 

OM 

O.2911 

00 

129.07 

-  47.42 

77-69 

00 

0.1814 

OQ 

54.00 

47.42 

81.42 

OQ 

O.OO83 

£ 

5108.48 

0  Li -99. 20 

4400.74 

L 

2.7449 

rj 

z 

9.965  X 

L0U  Lb -in*  (Applied  al  C.G.  ul*  A/ 

•a.  259-04) 

°Wln& 

Horlzuntal 

Stabilizer:  (Weight:;  D1 : 

itrlbutlon  4/17/62)  * 

Pane  L 

1  -Lb. In1 

uU 

Point  Wt.-Lb. 

H.S.10 

29.9 

6922 . 0 

H . S . 20 

21.5 

6024.0 

f» 

J* 

91.2 

14940.0 

Vertical.  Stub!  1 1  zer : 

(Welj'bt:;  DlatrlliuLUui  4/l8/<>2) 

* 

Component 

Basic 

v.:;.  h.s. 

Total  V.H. 

Basic  V.S. 

Pane  1 
Point 

WU- 

Lb.  Wt.-Lb. 

Wt.-Lb. 

1 

c 

-Lb-lri 

:u 

v.ii.  io 

.57.9 

57.9 

18974.0 

V.S.  PO 

PO. 

3 

20.5 

9456.0 

V.S. 

0 

01  .P 

91.2 

0 

>; 

‘18. 

01.  P 

109.4 

u 

'4410.0 

*  Ku^eli^'e  -  I /* '  Air[iltiiu  j  H»»ri  zontul  and  Vertical  stabilizer  -  Per  Side 
+  Dintrlbut ion  ui’  Lumped  Puurlof'e  Mann  at  Pane  1  PoLnt  01 
/  'lit  Applied  to  PunoLi*/*  fur  Projier  A  /(J  Yaw  inertiu  Simulation 


TABLE  6 


O  AILERON 

MASS  AND  INERTIA  CONFIGURATIONS 


Mass  Condition  1  -  Nominal:  (Weights  Distribution  4/18/62) 


Component 

Basic  Aileron  (Unbalanced) 

Maas  Balance’ 

Totul  Aileron  (Balanced) 

Panel  Point 

Wt.~  Lb. 

x~m. 

I|L~Lb.-ln.2 

Wt.~Lb. 

X~  In . 

I|L  Lb. -In. 2 

Wt.~Lb. 

X— In. 

1 1|  ~Lb.-ln.2 

55 

4.70 

4.14 

274.8 

a. 50 

-6.7 

160.0 

8.20 

-0.49 

433.8 

62 

2.14 

6.19 

190.5 

2.24 

-5.9 

78.0 

4.38 

0.01 

268.5 

70 

1.92 

5.48 

149.9 

2.22 

-5.3 

62.4 

4.14 

-0.30 

212.3 

76 

2.16 

5.18 

148.2 

2.38 

-4.7 

52.6 

4.54 

0 

200.8 

85 

2.48 

3.46 

181.2  1 

0.74 

-13.3 

130.9 

3.22 

-0.39 

312.1 

£  13.39  11.08  24.48 


£  WX  -  03.0  Lb. -In.  £WX=-69.5  £WX--6.5 


Mass  Condition  2:  (Balance  at  Outb'd.  Panel  Point) 


Component 

Mass  Balance’ 

Total  Aileron  (Balanced) 

Punel  Point 

W  -  Lb. 

X~in. 

WX—Lb.-In. 

I|L~Lb.-In.2 

W  -  Lb. 

X~ln. 

WX-Lb.-In. 

l|L~Lb.-ln.2 

55 

0 

0 

0 

0 

4.70 

4.14 

19.4 

274.8 

62 

0 

0 

0 

0 

2.14 

6.19 

13.2 

190.6 

70 

0 

0 

0 

0 

1.92 

5.48 

10.5 

149.9 

76 

0 

0 

0 

0 

2.16 

5.18 

11.2 

148.2 

85 

0.74 

-13,3 

-9.84 

130.9 

3.22 

-0.39 

-1.2 

312.1 

£ 

0.74 

-9.84 

14.14 

53.1 

Mass  Condition  3:  (Balance  at  Outb'd.  Panel  Point) 


Component 

Mass  Balance’ 

Total  Aileron  (Balanced) 

Panel  Point 

W  -  Lb. 

3(~In. 

WX~Lb.-ln. 

I,L~Lb.-In.2 

W  -  Lb. 

Y~In. 

WX~Lb.-In. 

I|L~U>.-In.2 

56 

0 

0 

0 

0 

4.70 

4.14 

19.4 

274.8 

62 

0 

0 

0 

0 

2.14 

6.19 

13.2 

190.5 

70 

0 

0 

0 

0 

1.92 

6.48 

10.6 

149.9 

76 

0 

0 

0 

0 

2.16 

6.18 

11.2 

148.2 

86 

6.29 

-10.0 

-62.9 

629.0 

8.77 

-6.19 

-64.2 

810.2 

£ 

6.29 

-62.9 

19.69 

0.10 

’Treated  as  a  Point  Mass 


TABLE  6  (CONTINUED) 


AILERON 

MASS  AND  INERTIA  CONFIGURATIONS 


Mass  Condition  4:  (Balance  at  Outb'd.  Panel  Point) 


Component 

Maas  Balance* 

Total  Aileron  (Balanced) 

Panel  Point 

Wt.~Lb. 

X~In. 

WX~Lb.-in. 

I^Lb.-In.2 

Wt.~Lb. 

X~ln. 

WX~Lb.-ln. 

1  „~Lb.-ln  2 
ft 

66 

0 

0 

0 

0 

4.70 

4.14 

19.4 

274.8 

62 

0 

0 

0 

0 

2.14 

6.19 

13.2 

190.6 

70 

0 

0 

0 

0 

1.92 

6.48 

10.5 

149.9 

76 

0 

0 

0 

0 

2.16 

6.18 

11.2 

146.2 

85 

9.62 

-10.0 

-96.2 

962.0 

12.10 

-7.24 

-87.6 

1143.2 

E 

9.62 

-96.2 

23.02 

-33.2 

Maes  Condition  &:  (Balance  at  Inb'd.  &  Outb'd,  Panel  Points) 


Component 

Muss  Balance* 

Total  Aileron  (Balanced) 

Punel  Point 

Wt.~Lb. 

Tr*  In. 

WX~Lb.-In. 

Wt.~U>. 

X~In. 

WX-Lb.ln. 

l^-Lb.-ln2 

65 

12.38 

-7.0 

-86.7 

606.6 

17.08 

-3.93 

-67.1 

881.4 

62 

0 

0 

0 

0 

2.14 

6.  19 

13.2 

190.6 

70 

0 

0 

0 

0 

1.92 

5.48 

10.5 

149.9 

76 

0 

0 

0 

0 

2.16 

5. 18 

11.2 

148.2 

86 

2.40 

-4.1 

-9.8 

40.3 

4.88 

-0.26 

-1.3 

221.6 

E 

14.78 

-96.5 

28.18 

-33.5 

Maes  Condition  6:  (Balance  at  Inb'd.  &  Outb'd.  Punel  Points) 


I 

Component 

Muse  Balance* 

Total  Aileron  (Ualunced) 

Punel  Point 

Wt.~U). 

X~In. 

WX~Lb.-ln. 

l1(~Lb.-ln.2 

Wt.~  Lb. 

X~Ln. 

w5<~I.U-la 

l|L~Lb.-!n.2 

65 

NO'I 

'  AVAILABLE 

62 

0 

0 

0 

0 

2.14 

6. 19 

13.2 

190.5 

70 

0 

0 

0 

0 

1.92 

5.48 

10.5 

149.9 

76 

0 

0 

0 

0 

2.  16 

6.  18 

11.2 

148.2 

85 

NOT  AVAILABLE  | 

♦Treated  uu  u  Point  Mass 


TABLE  6  (CONTINUED) 


AILERON 

MASS  AND  INERTIA  CONFIGURATIONS 


Mass  Condition  7:  (Balance  at  Inb'd.  &  Outh'd.  Panel  Points) 


Component 

Mass  Balance4 

Total  Aileron  (Hulanccii) 

Panel  Point 

Wt.  ~  Lb. 

WX~Lb.-ln. 

l,L~Lb.-ln.2 

Wt.~ib. 

X~ln. 

WX~Lb.-In. 

55 

8.53 

-7.0 

-59.71 

418.0 

13.23 

-3.04 

-40.2 

692.8 

62 

0 

0 

0 

0 

2.14 

6.19 

13.2 

190.5 

70 

0 

0 

0 

0 

1.92 

5.48 

10.5 

149.9 

76 

0 

0 

0 

0 

2.  16 

5.  18 

11.2 

148.2 

65 

2.4 

-4.1 

-9.84 

40.3 

4.86 

-0.25 

-1.2 

221.5 

£ 

10.93 

-69.5 

24.33 

-6.5 

Mass  Condition  8:  (Balance  at  Outb'd.  Panel  Point) 


Component 

Maas  Balance4 

Total  Aileron  (Balanced) 

Panel  Point 

Wt.  -  Lb. 

X~ln. 

WX-Lb.-In. 

Wt.~li>. 

Xr-ln. 

WX-Lb.-In. 

I„~L»>.-ln.2 

55 

0 

0 

0 

0 

4.70 

4. 14 

19.4 

274.8 

62 

0 

0 

0 

0 

2.14 

6.19 

13.2 

190.6 

70 

0 

0 

0 

0 

1.92 

5.48 

10  a 

149.9 

76 

0 

0 

0 

0 

2.16 

5.16 

U.2 

148.2 

85 

7.0 

-10. 0 

-70.0 

700.0 

9.48 

-6.47 

-61.3 

661.2 

£ 

7.0 

-70.0 

20.4 

-7.0 

Maas  Condition  9:  (Balance  at  Outb'd.  Panel  Point) 


Component 

Muss  Balance4 

Totul  Aileron  (Balanced) 

Panel  Point 

Wt.  -  Lb. 

X~ln. 

VVX~ll>.-ln. 

I  Lb. -In. 2 

Wt.~Lb. 

WX~LB.-In. 

I„;-Lb.-ln.2 

55 

0 

0 

0 

0 

4.70 

4. 14 

19.4 

274.  e 

62 

0 

0 

0 

0 

2.14 

6. 19 

13.2 

190.6 

70 

0 

0 

0 

0 

1.92 

5.48 

10.5 

149.9 

76 

0 

0 

0 

0 

2. 16 

6.18 

11.2 

148.2 

85 

5.0 

-10. 0 

-50.0 

500.0 

7.48 

-5.53 

-41.3 

661.2 

£ 

5.0 

-50.0  ' 

18.39 

13.0 

Note;  Unbalance  of  aileron  measured  with  respect  to  hinge-line,  positive  aft 


♦Treated  as  a  Point  Mass 


TABLE  7 


MODAL  DATA 


Stiffness  Configuration  OA 
Mass  Condition  9 
No  Gyroscopic  Forces 

Symmetric : 


Mode 

Frequency 

ops 

Dumpli^ 

Point  Driven 

Predominant  Component 

i 

6.82 

O.OltB 

LOA 

Fuse  luge 

2 

8.^0 

o.oi-o 

90K 

Wing 

5 

12.00 

0.09  5 

10A 

Fuselage 

4 

12.68 

0.066 

Ail.  By  (m)* 

Aile run -Fuselage 

6 

16.66 

0.056 

90K 

Wing 

6 

20.66 

0.022 

yoK 

Wing- Fuse  lage 

7 

25.96 

0.026 

yoK 

Wifig-Hurlz.  Stub!  Lizer 

8 

26.66 

0.021. 

H.S.  ^0 

Horizontal  Stabi 1 izer 

-llynum-trlo 

: 

Frequency 

Dumping 

Mode.- 

cps 

8 

Point  Driven 

Predominant  Component 

1 

8.  Bit 

o.otj7 

H.S.  r‘0 

Empennage 

') 

t 

10.6't 

0.0)0 

Ail.  By  (m)* 

Ai leron-Empennage 

5 

1  1  .t)6 

o.o  V* 

yoK 

Wing -Empennage 

It 

lit. Ol- 

u.aA 

H.S.  lJG 

Empennage 

6 

lO.  66 

0.007 

9 OK  (M)* 

Wing 

6 

20. ‘1 1. 

0.019 

OA 

Fuse*  1  age 

Y 

21.76 

o.ooB 

90K 

Wing 

8 

50. 10 

0.01') 

V.D.  50 

Fuse  luge -Empennage 

9 

Itl-.Ou 

o.o  ;*5 

6  OK 

Wing 

*  Moment 


TABLE  8 

MODAL  DATA 


Cti  lTnoiLi  Coni’i/'urntlon  9P 
Mar,:;  Condition  8 
No  Gyrotnruplc  Foivos 

Symmetric : 


Fivqiioin’y  Damping 

Modo  cpn  ^ _ Point  privon 


1 

O.oM 

o.o'o 

iOQ 

P 

M .  00 

O.OV.' 

90K 

I 

J 

11.44 

0  .uh‘j 

Ai  1 .  89  (M) 

It 

10.11 

O.ult  9 

10A 

* ) 

l‘j-w 

0.0)M 

9UK 

0 

.,'0.09 

o.o:^ 

90K 

7 

S-'.'Jlt 

o.o. '9 

9UK 

M 

U't.l  to 

0.0. ’4 

H.G.  PU 

Ant L-Symmotr i e ; 


k 

Frequency 

op:; 

Damp!  ig' 

8 

Point.  Driv<* 

1 

8.(9* 

U.u'jU 

HPk  PO 

P 

v.yo 

o.uoy 

Ail.  Mo  (M) 

i 

> 

i  1  .UP 

0.0';.' 

9UK 

it 

|4.V 

o.o.p. 

H.G.  PU 

i 

j 

10. MB 

0.0' 0 

1*00 

i, 

0.01H 

OA 

7 

.■o.M't 

0.010 

90K 

b 

.'V.Mo 

O.OlH 

OQ 

9 

111.  1  ‘ , 

0.0/'; 

UOK 

Pxvdominan t  Component 

Kutj*;  1  Uf'u 
W1  rg' 

Ai  ioron-Kunt*  la^o 
p\i : :  •  *  I  U*'t: 

WLn^ 

Winf'-Ku:J«*  1  u//o 
Wing-lloriz.  GtatiLl  iz*:r 
Horizontal  Gtalii  i  L/.or 


Prodorninmt  (Jumponun  t 

EmpniLiiut'c 

Ai  lnron-Kin£joniiUf'n 

Wintf-Emj*.-ntUi*n* 

Kmprnna^/o 

Win/1: 

Ku:>t;  1  u^r 
Win/1; 

Ku:;«  *  ]  afV?-Krnp«*nnu//t 1 
W  i  ntf 


*  Moment. 


TABLE  9 


MODAL  DATA 


Stiffness  Configuration  9A 
Muss  Condition  9 
Gyroscopic  Forces 

Symmetric : 


Mode 

Frequency 

cps 

Dumping, 

B 

Point  Driven 

Predominant  Component 

i 

0.65 

0.044 

L0A 

Fuse  luge 

» 

8.4  j 

0.041 

90K 

Wing 

5 

12.06 

0.030 

10A 

Fuse  luge 

4 

12.40 

0.060 

ALi.  89  (M)* 

Ai  leron-Fuselage 

9 

IO.34 

0.037 

90K 

Wing 

6 

20.40 

0.0. '9 

90K 

Wing-Fuselage 

7 

:*4.19 

0.026 

90K 

Wing-Horiz .  Stabilizer 

8 

29.66 

0.001 

li.n.  90 

Horizontal  Stabilizer 

1 

l  § 

r- 

Ki*<  •'iu<-n>,y 

Dcuti|i  i 

_  ...  a  - 

Point  Driven 

Predum  L  nan  t  Corn  pone  n  t 

l 

« .  Hi* 

0.0;  *8 

H.S.  90 

Empennage 

9 

10.  V 

0.099 

Ail  .  89  (M)» 

Ai  le ron -Empennage 

/ 

> 

i  i  .(>9 

0. 0.-50 

90K 

Wing -Empennage 

4 

1**  .00 

0.099 

H.S.  90 

Empennage 

\> 

10.  59 

0.0:  *9 

1*0G 

Wing 

(, 

90. 4o 

0.019 

OA 

Fuse  luge 

7 

pp .  00 

0.099 

9  OK 

Wing 

8 

w.09 

0.019 

OQ 

Fusel  age -Krrqje  linage 

9 

if,’.  00 

0.095 

90K 

Wing 

TABLE  10 


MODAL  DATA 


Stillness  Configuration  1>P 
Mass  Condition  8 
Gyroscopic  Forces 

Symmetric : 


Mode 

Frequency 

cps 

Dump lag 

JL  _ 

Point  Driven 

i 

6.05 

0.049 

10Q 

2 

8.00 

0.058 

90K 

3 

11.45 

0.047 

Ail.  89  (M)* 

it 

12.12 

0.049 

IGA 

3 

0-72 

0.039 

90K 

6 

20.12 

0.025 

90K 

7 

25.20 

0.026 

90K 

8 

29.»*9 

0.020 

H.S.  20 

Anti-Symmetric 

Mode 

Frequency 

cps 

Dumping 

6 

Point  Driven 

1 

8.64 

0.030 

H.S.  20 

2 

9.71 

0.097 

All.  89  (M)* 

3 

11.05 

0.032 

90K 

4 

14.90 

0.028 

H.S.  20 

19.72 

0.030 

40G 

6 

20.28 

0.019 

0A 

7 

20.96 

0.021 

90K 

8 

29.90 

0.017 

0Q 

9 

41.20 

0.023 

90K 

*  Moment 


Predominant  Component 

Fuselage 

Wing 

Ailerori-F  ;olage 
Fuse  log* 

Wing 

Wing-Fuse  lage 
Wing-Hori z.  Stabilizer 
Horizontal  Stabilizer 


Predominant  Component 

Empennage 

Ail  eron-Empennage 

Wing-Empennage 

Empennage 

Wing 

Fuselage; 

Wing 

Fuael  age -Empennage 
Wing 
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